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Abstract 


Multiple  use  productivity  is  described  for  watershed  lands  in  the 
ponderosa  pine  type  of  Arizona,  with  special  emphasis  on  the  Beaver 
Creek  Pilot  Watershed  near  Flagstaff.  Yields  of  timber,  herbage,  and 
water  under  past  management  are  reported,  along  with  information 
on  wildlife  values,  esthetics,  flood  and  sedimentation  hazards,  and 
water  quality.  Changes  in  productivity  and  environmental  quality 
are  then  described  following  five  experimental  land  treatments  in- 
cluding livestock  grazing  and  various  levels  of  forest  thinning  and 
clearing.  Preliminary  analytical  procedures  for  predicting  treat- 
ment responses  and  costs  allow  the  user  to  estimate  the  tradeoffs  in 
production  and  environmental  quality.  Some  further  research 
needed  for  bringing  these  and  other  response  models  to  operational 
capability  is  described. 
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Introduction 

This  report  documents  the  background  of 
water  yield  and  multiple  use  research  in  the 
ponderosa  pine  type  in  Arizona,  and  sets  forth 
the  objectives  of  the  program  initiated  on  the 
Beaver  Creek  Pilot  Watershed  on  the  Coconino 
National  Forest  near  Flagstaff  to  carry  out  the 
research.  It  characterizes  the  Beaver  Creek 
area  in  physical  terms,  and  describes  its  multi- 
ple use  production  under  present  management. 
This  section  is  followed  by  an  updated  presenta- 
tion of  results  of  the  various  treatments  being 
pilot  tested.  Models  are  described  that  eventual- 
ly will  be  incorporated  into  a  multiple  use 
management  planning  system  of  broad  appli- 
cability. Interrelationships  of  the  various 
products  are  discussed  and  implications  are 
pointed  out  indicating  the  best  general  forest 
density  levels  to  achieve  various  management 
objectives.  Finally,  needed  research  is 
described  to  provide  better  management  tools. 

Background 

In  the  mid-1950's,  water-user  and  ranching 
groups  in  the  Salt- Verde  Basin  called  attention 
to  the  possibility  of  treating  watersheds  to  in- 
crease irrigation  water  yields  and  provide  more 
grass  for  grazing.  In  response  to  these  con- 
cerns, the  Chief  of  the  Forest  Service  approved 
in  August  1956  the  "Arizona  Water  Management 
Project,  Beaver  Creek."  The  objectives  of  the 
project  were  to  observe  the  effects  of  vegeta- 
tive changes  on  water  yield,  soil,  forage,  fish, 
wildlife,  and  recreation,  and  to  determine 
whether  these  changes  increased  the  risk  from 
fire,  insects,  and  disease.  The  work  in  1956-57 
was  to  be  an  action  watershed  improvement 
program  involving  clearing  of  4,000  acres  of 
juniper  and  silvicultural  thinning  of  1,200  acres 
of  ponderosa  pine. 

During  this  period,  a  staff  of  watershed  man- 
agement experts  studied  the  economic  feasibil- 
ity of  obtaining  additional  water  for  irrigation 
from  the  watersheds  of  Arizona.  The  report  re- 
sulting from  this  study,  known  as  the  "Barr  Re- 


port," was  released  in  December  1956  (Barr 
1956). 

The  Arizona  Water  Resources  Committee  in 
February  1957  recommended  a  program  of  ac- 
tion and  research,  the  purpose  of  which  was  to 
find  out  if  modification  of  vegetation  could  pro- 
duce more  water  on  the  Salt-Verde  Watershed. 
Specific  recommendations  were  made  for  proj- 
ects in  spruce-fir,  ponderosa  pine,  juniper, 
chaparral,  and  streambed  vegetation  (Brown 
1970). 

The  Beaver  Creek  watershed  was  chosen  to 
represent  the  juniper  and  pine  types.  In  the 
period  1957-62,  stream  gages  were  built  on  18 
small  watersheds  of  66  to  2,036  acres:  three  in 
the  Utah  juniper  type,  three  in  the  alligator 
juniper  type,  and  12  in  the  ponderosa  pine  type. 
Stream  gages  were  also  built  on  two  large 
watersheds,  Woods  and  Bar  M,  to  provide  an 
opportunity  for  testing  composite  treatments 
developed  on  the  small  watersheds.  Sediment- 
measuring  installations  of  new  design  were 
constructed  on  nine  of  the  watersheds.  Also,  a 
network  of  precipitation  gages  was  installed. 
Timber,  range,  and  wildlife  were  inventoried, 
the  latter  by  the  Arizona  Game  and  Fish  De- 
partment. 

Objectives  of  the  Beaver  Creek 
Watershed  Evaluation  Program 

The  general  objective  of  the  project  as  as- 
signed in  the  1960's  was  to  evaluate  land  man- 
agement measures  designed  to  increase  water 
yields.  While  water  yield  was  to  receive  major 
consideration,  changes  in  forage  and  timber 
production,  wildlife  populations,  recreational 
value,  and  erosion  and  sediment  movement 
were  also  to  be  evaluated.  Costs  and  benefits  of 
the  various  treatments  were  to  be  studied  so 
that  the  total  project  could  be  evaluated.  These 
objectives,  maintained  until  1971,  resulted  in 
six  pilot  treatment  tests  on  pine  watersheds  and 
three  in  the  juniper. 

In  1971  an  Interdisciplinary  Team  assigned 
by  the  Chief  of  the  Forest  Service  reviewed  the 
Project  and  recommended  the  objectives  be 
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changed  to  read  as  follows: 

1.  To  develop  multiple  use  production  data  for 
alternative  land  management  practices,  and 
determine  their  ecological  and  environmen- 
tal impacts. 

2.  To  develop,  test,  and  refine  physical  product 
models  that  will  estimate  responses  to  man- 
agement activities  for  use  in  decisionmaking 
in  the  Southwest. 

3.  To  develop  and  test  an  economic  evaluation 
and  planning  model. 

4.  To  pilot-test  and  demonstrate  on  an  opera- 
tional scale  the  physical  product  and 
evaluation-planning  models  and  other  man- 
agement procedures  thus  developed. 

5.  To  develop  from  the  physical  and  economic 
models  a  series  of  management  models  with 
broader  geographic,  economic,  and  en- 
vironmental application  potential. 

6.  To  publish  a  report  explaining  the  physical 
responses  to  land  management  treatments, 
procedures  recommended  for  planning 
management  programs,  and  how  a  balanced 
land  management  program  in  the  Salt- Verde 
Basin  might  affect  resource  outputs  and  en- 
vironmental quality. 

This  report  will  document  the  progress  that 
has  been  made  toward  meeting  these  objectives. 
In  brief,  it  will  show  that  the  objectives  of  the 
sixties  and  the  first  objective  of  the  seventies 
have  largely  been  met.  The  modeling  and  large- 
scale  pilot  testing  objectives  of  the  seventies 
will  be  discussed,  although  only  limited  results 
can  be  presented. 

Physical  Characteristics  and  Productivity 
of  the  Ponderosa  Pine  Type 

Ponderosa  pine  forests  occupy  about 
4,282,000  acres  in  Arizona  and  4,681,000  acres  in 
New  Mexico  (Spencer  1966;  Choate  1966).  In  the 
Salt- Verde  Basin  in  Arizona  they  occupy  more 
than  1,650,000  acres  or  20  percent  of  the 
watershed  (fig.  1);  they  produce  nearly  one-half 
of  the  total  runoff  of  the  Salt  River  (Barr  1956). 
The  Beaver  Creek  watershed,  to  be  charac- 
terized in  this  section,  is  about  one-third  pon- 
derosa pine  (fig.  2).  Total  area  of  Beaver  Creek 
is  about  275,000  acres. 

Physiography  of  the  Pine  Type 
on  Beaver  Creek 

A  high  plateau,  sloping  mesas  and  breaks, 
steep  canyons,  and  valleys  characterize  the  to- 
pography on  Beaver  Creek.  Elevations  in  the 
pine  type  vary  from  about  6,800  feet  to  8,000 
feet.  The  bedrock  underlying  the  area  consists 
of  igneous  rocks  of  volcanic  origin;  below  them 
are  sedimentary  rocks  of  Kaibab,  Coconino,  and 
Supai  formations. 


Figure  1 . — The  ponderosa  pine  type  in  Arizona. 


Geology  of  the  watershed  was  studied  by 
Beus,  Rush,  and  Smouse,^  and  Scholtz  (1968). 
Their  inventories  included,  in  addition  to  rock 
types,  identification  of  geologic  features  such 
as  faults  and  fractures  that  can  be  quantified 
and  used  in  water  prediction  models. 

The  soils,  developed  on  basalt  and  cinders,  are 
mostly  silty  clays  and  silty  clay  loams  less  than 
2-1/2  feet  deep.  They  are  described  in  detail  by 
Anderson,  Williams  and  Creeze  (1960)  and  Wil- 
liams and  Anderson  (1967).  To  better  predict  the 
hydrologic  properties  of  the  soils,  the  relation- 
ships between  physical  and  hydrologic  charac- 
teristics of  soils  within  the  various  soil  man- 
agement areas  were  studied  in  the  laboratory  by 
Rector  (1969),  Ryan  (1969),  and  Thames.^ 

Hydrologic  Characteristics  of  the 
Pine  Type  on  Beaver  Creek 

This  section  briefly  describes  the  hydrologic 
characteristics  of  the  southwest  ponderosa  pine 
type,  and  in  particular  the  Beaver  Creek 

"Unpublished  reports  available  at  the  Rocky  Mt.  For. 
and  Range  Exp.  Stn.,  Flagstaff,  Ariz. 

•Thames,  John  L.  A  study  to  determine  the  hydrologic 
and  physical  properties  of  some  Beaver  Creek  soils.  Final 
report  of  a  cooperative  study  between  the  Univ.  of  Ariz, 
and  the  Rocky  Mt.  For.  and  Range  Exp.  Stn.,  331  p.  Flags- 
taff, Ariz.  1969. 
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Figure  2. — The  Beaver  Creek  Pilot  Watershed. 


watershed.  It  begins  with  a  general  description 
of  the  temperature  and  precipitation  regimes, 
and  then  follows  precipitation  as  it  falls  on  a 
watershed  to  its  eventual  disposition  as  runoff, 
evaporation,  transpiration,  or  subsurface  flow. 

Temperature. — The  Beaver  Creek  watershed 
lies  in  the  center  of  Arizona  in  what  is  called  the 
plateau  climatic  region,  which  includes  the 
Kaibab,  Coconino,  and  MogoUon  Plateaus  and 
the  San  Francisco  and  White  Mountains  (Green 
and  Sellers  1964).  January  temperature  in  the 
region  averages  26°  F  and  the  July  tempera- 
tures 66"  F  (Schubert  1974). 

Climatic  records  of  different  kinds  are  pres- 
ently being  obtained  at  50  locations  in  the  pine 
type  on  Beaver  Creek.  At  ten  of  the  locations 
precipitation  is  recorded  continuously.  Temp- 
erature and  humidity  are  recorded  at  three  lo- 


cations. Average  annual  temperature  in  the  pine 
type  has  been  45°  F.  Monthly  averages  have 
been  as  follows: 


Maximum  Minimum 

 (Degrees  F.) 


Mean 


October 

64 

30 

47 

November 

52 

23 

37 

December 

44 

15 

29 

January 

44 

14 

29 

February 

45 

16 

30 

March 

49 

19 

34 

April 

56 

24 

40 

May 

67 

31 

49 

June 

76 

37 

57 

July 

83 

48 

66 

August 

78 

47 

63 

September 

73 

41 

57 
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These  data  are  based  on  an  average  of  12 
years  of  record  from  weather  stations  on 
Watersheds  8,  17,  and  Bar  M. 

Precipitation. — The  ponderosa  pine  type  in 
Arizona  receives  between  15  and  25  inches  of 
precipitation  annually  (Schubert  1974).  On 
Beaver  Creek  the  average  annual  precipitation 
has  ranged  from  18  to  35  inches  with  an  average 
of  25  inches  on  six  representative  watersheds. 
These  watershed  averages,  based  on  13  years  of 
record,  were  determined  by  the  Theissen 
method  using  20  precipitation  gages.  Precipita- 
tion distribution  throughout  the  water  year  is 
illustrated  by  these  monthly  averages  for 
Watersheds  7,  8,  10,  13,  15,  and  18: 


Inches 

October 

1.08 

November 

2.40 

December 

3.65 

January 

2.10 

February 

2.25 

March 

2.59 

April 

1.83 

May 

.65 

June 

.44 

July 

2.34 

August 

3.22 

September 

2.41 

Average  annual 

24.96±0.61 

There  are  two  major  precipitation  seasons  in 
the  plateau  region.  The  most  important  is  dur- 
ing the  winter  —  October  through  April.  Sixty- 
four  percent  of  the  annual  precipitation  falls 
during  these  7  months  on  Beaver  Creek.  The 
second  is  during  the  summer,  particularly  July 
through  September,  during  which  32  percent  of 
the  annual  precipitation  falls. 

Predictable  atmospheric  circulation  patterns 
carry  large  quantities  of  moisture  over  the  State 
and  encourage  precipitation.  Moisture  for  most 
summer  precipitation  comes  from  the  Gulf  of 
Mexico.  Occasional  record  summer  rains  in  the 
past  century  have  been  associated  with  move- 
ments of  tropical  air  into  the  State  from  the  Gulf 
of  California  and  the  Pacific  Ocean.  Moisture 
for  winter  precipitation  comes  from  the  Pacific 
Ocean. 

Average  annual  snowfall  for  15  locations  in 
the  ponderosa  pine  forests  of  Arizona  ranges 
from  a  low  of  12  inches  to  a  high  of  94,  with  an 
annual  mean  of  46  inches  (Schubert  1974). 
Snowfall  is  not  consistent  from  year  to  year  in 
this  area.  Fort  Valley  averages  91  inches,  how- 
ever in  13  out  of  60  winter  seasons  less  than  60 
inches  of  snow  fell.  Snowfall  during  34  percent 
of  these  winters  averaged  less  than  1  foot  per 
month  for  four  consecutive  months. 


The  snow  regime  as  described  by  the  SCS 
snow  courses  at  Casner  Park,  Happy  Jack,  and 
Munds  Park  provide  a  usable  estimate  of  the 
snow  conditions  on  Beaver  Creek.  The  15-year 
average  water  equivalents  at  these  locations  on 
March  1  are  3.2,  4.4,  and  2.7  inches,  respec- 
tively, or  3.4  inches  overall.  This  compares  to 
3.6  inches  at  Heber  and  Workman  Creek, 
Arizona.  The  snow  regimes  on  Beaver  Creek  are 
characteristic  of  locations  along  the  Mogollon 
Rim;  however,  they  are  higher  than  pine  areas 
away  from  the  Rim  such  as  Fort  Valley,  which 
has  an  average  of  2.1  inches. 

Snow  regimes  on  Beaver  Creek  are  quite  var- 
iable and  result  in  different  runoff  patterns  and 
amounts.  Most  runoff  is  produced  when  there  is 
a  continuous  snowpack  from  November 
through  April.  In  these  years  there  is  often  a 
continuous  small  flow  of  water  from  the  pine 
watersheds,  followed  by  a  large  runoff  period  in 
March  or  April.  Often,  major  runoff  occurs 
when  rain  falls  on  a  snowpack.  In  Water  Year 
1966  there  was  almost  continuous  flow  all 
winter  as  a  result  of  intermittent  storms  which 
dropped  mixed  precipitation  —  both  rain  and 
snow  —  on  a  shallow  snowpack.  In  Water  Year 
1972  a  snowpack  started  to  accumulate  in 
November  and  built  up  till  late  in  December 
when  a  rainstorm  removed  all  of  the  snow,  pro- 
ducing a  large  runoff  event.  A  third  type  of  snow 
regime  is  characterized  by  intermittent  snow- 
fall followed  by  a  dry  period  in  which  the  snow- 
pack disappears.  These  typically  light  snowfall 
years  yield  little  if  any  runoff  since  most  of  the 
water  in  the  snow  evaporates  back  to  the  atmos- 
phere or  infiltrates  into  the  soil.  In  a  fourth  snow 
regime,  very  little  snowfall  occurs  throughout 
the  season.  During  these  winters  the  snowpack 
seldom  accumulates  over  a  few  inches,  and  re- 
sults in  little  if  any  runoff. 

Interception,  stemflow,  and  throughfall. — 

Studies  of  rainfall  interception  by  pole-sized 
ponderosa  pine  in  central  Arizona  have  shown 
that  11  to  25  percent  of  the  summer  rainfall  is 
intercepted  (U.S.  Forest  Service  1958).  Of  the 
net  rainfall  (throughfall  plus  stemflow),  71  to  88 
percent  reaches  the  soil  as  throughfall  and  1  to 
8  percent  is  transported  to  the  soil  as  stemflow. 
Although  little  is  known  about  interception  of 
snow  by  ponderosa  pine  on  Beaver  Creek, 
Kittredge  (1948)  reported  22  percent  intercep- 
tion by  mature  ponderosa  pine  in  Idaho  and  18 
percent  for  young  ponderosa  pine  in  Colorado 
(Kittredge  considered  interception  to  include 
stemflow).  Thus  ponderosa  pine  apparently 
intercepts  similar  percentages  of  snow  and 
rain. 

Evaporation. — Studies  of  evaporation  from 
forest  floor  litter  under  ponderosa  pine  in  cen- 


4 


tral  Arizona  show  that  between  7  and  25  percent 
of  the  total  summer  rain  is  evaporated  to  the 
atmosphere  (U.S.  Forest  Service  1959). 

Studies  of  evaporation  from  snow  under  pon- 
derosa  pine  in  Arizona  and  New  Mexico  show 
that  between  0.006  and  0.16  inch  of  water  per  day 
is  lost  to  evaporation  (U.S.  Forest  Service  I960). 
Thus  a  snow  cover  present  for  2  months  may 
lose  as  much  as  5  inches  of  water  by  evapora- 
tion. In  open  areas  within  the  forest,  this  amount 
would  probably  be  greater. 

Infiltration. — The  major  soil  series  in  the 
ponderosa  pine  area  on  Beaver  Creek  are  Broll- 
iar,  Siesta,  and  Sponseller.  These  soils  are  de- 
veloped on  parent  materials  of  basalt  and  vol- 
canic cinders.  Their  textures  range  from  mod- 
erately fine  to  fine,  and  they  characteristically 
have  a  slow  rate  of  infiltration  —  0.05  to  0.2  inch 
per  hour  when  thoroughly  wetted  (Williams  and 
Anderson  1967).  These  soils  also  have  a  layer 
that  impedes  downward  movement  of  water. 
These  characteristics  are  major  factors  in  pro- 
ducing the  surface  runoff  which  comprises 
most  of  the  water  yield  on  the  Beaver  Creek 
watersheds. 

Streamflow. — New  techniques  in  streamgag- 
ing  have  been  utilized  on  Beaver  Creek  to  accu- 
rately measure  the  broad  range  of  debris-laden 
discharges  carried  by  the  streams.  A  concrete 
trapezoidal  flume  was  designed  for  use  on  the  18 
small  watersheds.  The  design  was  a  modifica- 
tion of  a  Washington  State  College  flume, 
adapted  to  Beaver  Creek  conditions  (fig.  3).* 
Another  was  a  new  design  for  gaging  discharges 
of  several  thousand  ft^/s  (cubic  feet  per  second) 
on  Woods  Canyon  and  Bar  M  (fig.  4;  Brown 
1969).  H  flumes  have  been  constructed  on 
24  small  homogeneous  subwatersheds  where 
flows  are  usually  less  than  15  ft^/s  (fig.  5). 

The  watersheds  are  organized  into  an  ex- 
perimental design  to  provide  a  basis  for  statisti- 
cal evaluation  of  changes  in  streamflow.  Water- 
sheds within  each  vegetation  type  are  paired, 
with  some  scheduled  for  treatment  and  others 
retained  untreated  as  statistical  controls.  Fol- 
lowing pretreatment  calibration,  one  watershed 
in  each  pair  is  treated.  These  methods  are  de- 
scribed by  Kovner  (1956),  Kovner  and  Evans 
(1954),  and  Wilm  (1943). 

Average  annual  streamflow  from  all  un- 
treated pine  watersheds  on  Beaver  Creek  has 
varied  from  0.2  to  16.4  inches  with  an  overall 
average  for  the  type  of  5.3±0.5  inches  per  year 
(based  on  14  years  of  record,  through  Water 
Year  1973).  Ninety-three  percent  of  this  flow 
has  come  during  the  period  October  through 
April  as  a  result  of  snowmelt  or  winter  rains. 

''Robinson,  A.R.  Study  of  the  Beaver  Creek  measuring 
flumes.  Report  to  the  Rocky  Mt.  For.  and  Range  Exp.  Stn. 
14  p.  Flagstaff,  Ariz. 


Figure  3. — Concrete  trapezoidal  flume  of  the  type 
used  on  Watersheds  1-18  on  Beaver  Creek. 


Figure  4. — The  Woods  Canyon  streamgage. 


Figure  5. — Streamgaging  installation  used  on  sub- 
watersheds. 
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Occasional  summer  storms  produce  flood  flows 
of  short  duration.  Although  these  floods  do  not 
contribute  much  to  the  total  water  yield  they  are 
a  major  factor  in  erosion  and  sedimentation. 

Figure  6  illustrates  the  monthly  distribution 
of  streamflow  and  precipitation  (based  on  112 
station  years  of  record).  For  the  year  as  a  whole, 
20  percent  of  the  precipitation  becomes 
streamflow.  For  the  October- April  period, 
streamflow  is  28  percent  of  precipitation,  and 
for  May-September  it  is  3  percent. 


PRECIPITATION 


STREAMFLOW 

I     t  L 


OCT      NOV      DEC      JAN      FEB      MAR      APR      MAY       JUN      JUL      AUG  SEP 

Figure  6. — Distribution  of  precipitation  and  stream- 
flow  on  the  Beaver  Creel<  ponderosa  pine  water- 
sheds. 

Sediment. — Sediment  yield  data  have  been 
obtained  by  means  of  catchment  basins  and 
splitting  devices  installed  on  Watersheds  9,  10, 
12, 14, 16,  and  17.  These  installations  (fig.  7)  are 
described  in  detail  by  Brown  et  al.  (1970).  Sus- 
pended sediment  concentrations  have  also  been 
sampled  intensively  on  all  watersheds.  Approxi- 
mately 1,000  samples  have  been  analyzed,  and 
the  results  are  to  be  reported  in  a  later  publi- 
cation. 

Mean  annual  sediment  yields  from  untreated 
pine  watersheds  (26  station  years)  have  varied 
from  a  trace  (less  than  0.005  ton  per  acre)  to  6 
tons  per  acre,  with  an  average  of  0.6  ton  per  acre 
over  an  11-year  period  of  record.  This  average, 
however,  should  not  be  used  to  characterize 
total  yields  from  pine  areas  because  it  is  influ- 
enced heavily  by  the  100-year  storm  of  Septem- 
ber 5,  1970.  Excluding  this  extreme  event, 
annual  sediment  yields  have  averaged  only  0.02 
ton  per  acre. 

Except  for  the  September  1970  storm,  essen- 
tially all  sediment  has  been  produced  during 


winter  seasons.  Winter  sediment  yields  have 
ranged  from  a  trace  to  0.06  ton  per  acre  with  an 
average  of  0.02  ton  per  acre. 

Sediment  yields  from  individual  untreated 
watersheds  have  ranged  from  a  trace  to  11  tons 
per  acre  during  the  11-year  period  of  record. 
The  11  tons  per  acre  were  derived  from 
Watershed  10  during  the  September  1970  storm 
which  produced  a  peak  streamflow  discharge  of 
1,090  ft^/s/mi^  with  an  estimated  recurrence 
interval  of  over  200  years. ^  This  watershed  was 
within  or  near  the  center  of  maximum 
30-minute  precipitation  intensity  and  total 
storm  precipitation,  3.0  inches  per  hour  and  7.0 
inches,  respectively. 

The  only  peak  discharge  recorded  on  Beaver 
Creek  which  exceeds  the  amount  in  this  storm 
was  from  a  3-year-old  clearcut  on  Watershed  12 
which  had  a  discharge  of  1,200  ft^/s/mi^and  a 
sediment  yield  of  9  tons  per  acre.  Clearcutting 
in  this  case  appears  responsible  for  the  high 
peak  discharge  and  resulting  high  sediment 
yield. 

On  untreated  Watershed  10  the  high  stream- 
flow  was  apparently  due  to  the  high  precipitation 
intensities.  Contributing  further  to  the  high  sed- 
iment yield  was  the  apparent  accumulation  of 
sediment  available  in  its  channels;  during  its 
7-year  period  of  record  it  had  never  produced 
more  than  0.03  ton  per  acre.  Therefore,  poten- 
tial sediment  yields  of  at  least  11  tons  per  acre 
may  be  realized  from  individual  untreated  pine 
watersheds  with  characteristics  similar  to 
those  found  on  Beaver  Creek. 

^Based  on  estimates  made  by  the  U.S.  Geological  Sur- 
vey. 


Figure  7. — Installation  for  obtaining  sediment  yield 
data.  Coarser  sediments  are  caught  in  the  basin 
above  the  dam,  while  suspended  sediment  is  sam- 
pled continuously  by  a  series  of  three  splitters. 
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Water  Quality 

Quality  of  runoff  water  has  been  analyzed 
since  1969.  Only  occasional  samples  have  been 
collected  and  sampling  conditions  have  varied. 
Most  samples  were  of  runoff  from  a  melting 
snowpack.  An  exception  was  a  group  of  four 
samples  from  the  September  5,  1970,  storm. 
Samples  were  taken  at  the  flumes,  either  dipped 
or  with  a  DH-48  sampler. 

Average  values  for  key  water  quality  charac- 
teristics on  untreated  pine  watersheds  are 
summarized  as  follows: 

Ion  Concentration 

Mg/1 


TDS 

44 

Ca 

5.1 

Na 

1.9 

PO4 

.16 

NO3 

.16 

Fe 

.29 

Electrical  conductivity  ranged  from  38  to  65  ix 
mho/cm.  Sodium  was  low,  with  SAR  (sodium  ab- 
sorption ratio)  below  0.2  on  all  samples.  No  one 
watershed  produced  water  of  consistently  dif- 
ferent quality  than  any  other.  Water  from  the 
pine  watersheds  met  quality  standards  for 
drinking  water,  aquatic  life,  and  irrigation  (Na- 
tional Technical  Advisory  Committee,  FWPCA 
1968).  Only  iron  approached  the  unacceptable 
level  for  drinking  water. 

Timber 

A  comprehensive  description  of  the  timber 
stands  on  Beaver  Creek  is  essential  to  this  pres- 
entation, because  results  of  the  treatments  and 
predictions  produced  by  the  current  models  are 
most  applicable  to  these  kinds  of  forest  stands. 

The  forest  on  Beaver  Creek  is  composed  of  a 
mixture  of  ponderosa  pine  (Pinus  ponderosa 
Laws.),  Gambel  oak  {Quercus  gambelii  Nutt.), 
and  alligator  juniper  (Juniperus  deppeana 
Steud.)  with  a  combined  basal  area  of  115 
ft^/acre.  Pine  basal  area  averages  92  ft^  Gambel 
oak  18  ft^  and  alligator  juniper  5  ft^  Volumes  of 
typical  untreated  stands  average  2,050  ft'^/acre 
overall,  of  which  1,775  ft^  are  pine  (Myers  1973), 
230  ft^  Gambel  oak  (Gevorkiantz  and  Olsen 
1955),  and  45  ft^  juniper  (Howell  and  Lexen 
1939).  Average  board-foot  volume  of  ponderosa 
pine  is  4,070  fbm/acre  (FfoUiott  et  al.  1971)  of 
which  only  5  percent  are  in  grade  3  saw  logs  or 
better.  Average  growth  of  pine  stands  is  40  ft^  or 
75  fbm/acre/year. 

The  number  of  trees  per  acre  by  size  class 
(table  1)  is  important  because  it  characterizes 
the  typical  stand.  There  is  an  excess  of  trees  in 


the  2-  through  10-inch  and  26-inch  plus  classes, 
and  a  shortage  of  trees  in  the  14-  through  24-inch 
classes.  Average  site  index  (height  in  feet  at  100 
years)  is  60  (Meyer  1938). 

Table  1. --Average  number  of  trees  per  acre,  by 
species  and  diameter  class,  on 
Beaver  Creek  watersheds  S- 1 0  and 
12-18 

D  iameter    pQp_  Al 1  i- 

class       ,  Gambel     qator    Quaking  All 

,        derosa  ,         .     .  . 

(Inches        .  oak       juni-      aspen  species 

J  ,   ,    ,     p I ne  1 

a.  b.n.  J     ^  per  ^  

------  Exmber  ------ 

Repro- 


ct  ion 

174. 

67 

82. 

91 

4.95 

262.53 

i. 

3^*9.83 

/in 

Dp 

1  Q  I'i  n 
\y . / i  u 

Z.H 

n 1 U . Hp 
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90. 

88 

20. 

1  3 

4.98 

03 

116. 02 

6 

51. 

09 

12. 

66 

1  .05 

02 

64.82 

8 

36. 

66 

5. 

46 

.38 

42.50 

10 

23. 

66 

3. 

73 

.10 

27.49 

12 

12. 

39 

2. 

99 

.12 

15.50 

lif 

6. 

31 

1 

76 

.05 

8.12 

16 

2. 

60 

1  . 

03 

.12 

3.75 

18 

2. 

02 

74 

.09 

2.85 

20 

1  . 

32 

47 

.07 

1 .86 

22 

1 . 

21 

20 

.  OH 

1 .45 

Ik 

97 

11 

.09 

1.17 

26 

70 

06 

.07 

.83 

28 

52 

02 

.07 

.61 

30 

28 

01 

.05 

.34 

32 

13 

.02 

.15 

07 

.03 

.10 

36 

03 

.Oh 

.07 

38 

01 

.ok 

.05 

01 

.01 

.02 

kl 

.01 

.01 

kk 

.02 

.02 

if6 

.01 

.01 

i»8 

.01 

.01 

50 

.02 

.02 

)tal2 

580 

69 

90 

.02 

n  .11 

.29 

698.22 

■  Peroent^- 

83 

17 

12 

89 

3.9 

.Ok 

100.00 

Includes  all  Juniperus  species. 
Exclusive  of  reproduction. 


The  inventory  system  used  on  Beaver  Creek 
went  beyond  the  standard  estimates  of  volume, 
basal  area,  and  number  of  trees  per  acre,  to 
include  timber  quality  data  describing  basic 
stem  characteristics  and  defect  features  that 
influence  quantity  and  quality  for  most  primary 
products  (Barger  and  Ffolliott  1970).  Further 
refinements  allow  levels  of  stocking  (Ffolliott 
and  Worley  1973)  to  be  calculated,  such  as  the 
percent  of  the  area  stocked  to  50  f  t^  basal  area  in 
trees  18  inches  or  larger.  From  these  stocking 
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levels,  timber  stocking,  recreation  potential,  or 
forage  production  can  be  determined  (Ffolliott 
and  Worley  1965). 

The  number  of  trees  per  acre  and  site  index 
are  used  to  initialize  the  timber  simulation 
model  described  in  a  later  section.  Growth  in- 
formation will  be  used  to  validate  model  predic- 
tions. Basal  area  is  the  stand  density  parameter 
used  in  production  functions  because  of  its  ef- 
fects on  other  products  such  as  runoff,  sedi- 
ment, and  forage  yields. 

Range 

Total  herbage  production  in  the  pine  type  on 
Beaver  Creek  averages  198  pounds  per  acre  an- 
nually for  untreated  conditions,  consisting  of 
124  pounds  of  perennial  grasses,  58  pounds  of 
forbs  and  half-shrubs,  and  16  pounds  of  shrubs. 
Following  are  the  major  herbaceous  and 
shrubby  species  occurring  under  ponderosa 
pine  stands  (average  pretreatment  values  for 
Watersheds  8,  9, 10, 12, 13, 14,  and  17  based  on  2 
to  9  years  of  record): 

Grasses  and  grasslike  plants   Pounds  per  acre 

Bluegrass,  mutton 

(Poafendleriana)  23 
Dropseed,  black 

(Sporobolus  interruptus)  10 
Grama,  blue 

(B out eloua  gracilis)  23 
Sedge  (Carex  spp.)  6 
Squirreltail,  bottlebrush 

(Sitanion  hystrix)  46 
Others  16 

Forbs  and  half-shrubs 

Aster,  showy 

(Aster  commutatus)  3 
Fleabane,  spreading 

(Erigeron  divergens)  4 
Fleabane,  trailing 

(E.  flagellaris)  3 
Goldeneye,  showy 

(Viguiera  multiflora)  3 
Ragweed,  western 

{Ambrosia  psilostachya)  7 
Snakeweed,  broom 

(Giitierrezia  sarothrae)  3 
Wormwood 

(Artemisia  carruthii)  3 
Others  32 

Shrubs 

Locust,  New  Mexican 

(Robinia  neomexicana)  3 
Oak,  Gambel 

(Quercus  gambelii)  13 


Total  198 


The  range  inventories  on  Beaver  Creek  origi- 
nally involved  only  transect  measurements  of 
condition  and  trend,  and  plot  measurements  of 
production  and  utilization.  Such  measurements 
have  been  made  periodically  at  10  locations 
within  Watersheds  7  through  14.  At  each  loca- 
tion there  is  a  cluster  of  ten  9.6-ft2  herbage  pro- 
duction and  utilization  plots  and  two  100-foot 
Parker  3-step  transects.  Herbage  production  is 
determined  for  each  species  by  weight  estimate 
and  checked  by  clipping  and  weighing 
(Pechanec  and  Pickford  1937a).  Forage  utiliza- 
tion is  determined  by  the  ocular-estimate-by- 
plot  method  (Pechanec  and  Pickford  1937b). 
Ground  cover  and  plant  frequency  are  deter- 
mined by  the  method  of  Parker  (1954).  The 
Parker  transects  aid  in  determining  the  ecolog- 
ical impact  of  the  different  treatments  on  plant 
composition  and  soil  stability. 

These  conventional  methods  of  range  inven- 
tory, however,  do  not  lend  themselves  to  statis- 
tical sampling  of  large  diverse  areas.  Thus  it 
became  necessary  to  measure  herbage  produc- 
tion at  the  timber  inventory  plot  locations  and  to 
develop  re.'ationships  between  herbage  produc- 
tion and  more  easily  measured  site  characteris- 
tics of  tree  density  and  cutting  history,  soils, 
topography,  forest  floor  characteristics,  and 
composition  and  distribution  of  the  forage. 
Some  of  these  relationships  are  described  by 
Clary  (1964,  1969);  Clary  and  Ffolliott  (1966, 
1969);  Clary,  Ffolliott,  and  Jameson  (1968); 
Clary,  Ffolliott,  and  Zander  (1966);  and  Clary 
and  Pearson  (1969)  (see  later  section  on  model- 
ing). In  addition,  special  studies  have  been  in- 
itiated to  characterize  range  response  following 
unique  treatments  such  as  stripcutting  and 
burning. 

Wildlife 

Big  game  use  on  Beaver  Creek  is  measured  in 
cooperation  with  the  Arizona  Game  and  Fish 
Department.  Changes  in  the  deer  and  elk  use  are 
estimated  from  counts  of  fecal  droppings  on 
plots  superimposed  on  the  timber  inventory 
sample.  Forage  preferences  of  deer  are  deter- 
mined by  observing  the  plants  eaten  by  tame 
deer  under  different  treatment  conditions.  Pro- 
duction of  preferred  wildlife  forage  is  deter- 
mined by  combining  species  yields  and  prefer- 
ence information.  Observational  studies  have 
been  conducted  to  determine  the  response  of 
wild  deer  to  different  treatment  situations. 
Hunter  check  stations  and  traffic  counters  have 
yielded  information  on  the  social  and  economic 
significance  of  the  game  resource. 

Big  game  numbers  have  declined  drastically 
on  Beaver  Creek  throughout  the  period  of  the 
study,  as  they  have  in  Arizona  and  the  South- 
west as  a  whole.  Deer  and  elk  populations  calcu- 
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lated  from  Neff's  1959-60  data  (Neff  1972) 
were  9  and  3  per  section,  respectively.  During 
the  period  1965  to  date  deer  populations  on  un- 
treated watersheds  have  been  approximately 
0.9  per  section  and  elk  1.3  per  section. 

The  forage  plants  on  untreated  areas  that  are 
most  preferred  by  deer  include,  in  order  of  pref- 
erence (Neff  1974):  „  J. 

Percent  of  deer  diet 

Gambel  oak 

(Quercus  gambelii) 
Unidentified  immature 


grasses 
Ponderosa  pine 

(Pinus  ponderosa) 
Redroot  eriogonum 

(Eriogonum  racemosum) 
Geranium  (Geranium  spp.) 
Clover  (Trifolium  spp.) 
Milkvetch 

(Astragalus  recurvus) 
Mountainmahogany 

(Cercocarpus  spp.) 
Red-and-yellow  pea 

(Lotus  wrightii) 


38 

7 

5 

4 
4 
3 

2 

2 

2 

67 


Because  of  concern  over  environmental  is- 
sues, the  wildlife  program  on  Beaver  Creek  was 
expanded  in  1972  to  include  studies  of  birds, 
insects,  carnivores,  and  desert  cottontail, 
Abert's  squirrel,  and  other  forest  rodents.  The 
ultimate  goal  will  be  to  describe  the  interactions 
between  all  major  forest  fauna  and  other 
ecosystem  components  under  a  variety  of  man- 
agement conditions.  The  latter  research  is 
being  conducted  in  cooperation  with  Northern 
Arizona  University,  the  University  of  Arizona, 
and  the  Museum  of  Northern  Arizona. 

Esthetics 

Much  of  the  esthetic  quality  of  forest  lands 
can  be  attributed  to  the  physiographic,  climatic, 
and  vegetative  characteristics  of  the  area.  The 
ponderosa  pine  forests  of  Arizona  are  situated 
in  an  intermediate  position  between  the  State's 
widespread  desert  and  woodland  regions  and  its 
higher  (but  widely  separated)  montane  forests 
of  aspen  and  mixed  conifers.  Because  this  in- 
termediate zone  has  a  vigorous  climate  with  rel- 
atively large  seasonal  variation,  it  provides  a 
welcome  visual  relief  for  travelers  and  re- 
creationists  from  other  parts  of  the  State.  Its 
potential  as  an  area  of  high  scenic  value  was 
suggested  by  an  1858  military  report: 

It  is  the  most  beautiful  region  I  ever  re- 
member to  have  seen  in  any  part  of  the  world. 
A  vast  forest  of  gigantic  pine,  intersected  fre- 


quently with  open  glades,  sprinkled  all  over 
with  mountains,  meadows  and  wide  savan- 
nahs, and  covered  with  the  richest  grasses, 
was  traversed  by  our  party  for  many  days  (Lt. 
E.F.  Beale,  quoted  by  Cooper  1960). 

The  regional  geology  is  such  that  the  south- 
trending  formations  of  the  Colorado  Plateau  di- 
rect groundwater  to  feed  several  perennial 
streams  which  begin  in  an  otherwise  arid  envi- 
ronment below  the  Mogollon  Rim.  Above  the 
Rim  on  the  plateau,  a  number  of  well-defined 
but  in  some  cases  ephemeral  streams  begin, 
which  feed  either  the  Little  Colorado  River  or 
the  Salt- Verde  systems. 

The  Beaver  Creek  study  area  includes  water- 
sheds in  both  the  woodland  and  the  intermediate 
zones.  Although  the  pine  watersheds  are  above 
the  Rim  topographically,  they  nevertheless 
drain  southward  to  the  Salt-Verde  River  system 
and  all  of  them  contain  ephemeral  streams. 

Litton  (1968)  developed  the  concept  of 
characterizing  landscapes  or  segments  of  land- 
scapes by  means  of  definable  terms  such  as 
"panoramic,"  "feature,"  "enclosed,"  "focal," 
"canopied,"  "detail,"  and  "ephemeral."  Land- 
scapes can  be  also  described  as  either  "macro," 
extending  far  into  the  distance,  or  "micro," 
where  visibility  is  restricted  to  a  matter  of 
yards  (U.S.  Forest  Service  1973). 

Using  this  terminology  we  can  briefly  charac- 
terize the  landscape  of  the  ponderosa  pine  type 
on  Beaver  Creek,  while  considering  at  the  same 
time  the  possible  implications  of  several  forest 
management  practices. 

These  lands  generally  afford  few,  if  any, 
panoramic  landscapes.  They  are  too  flat,  too  cov- 
ered with  forest  trees  to  provide  vistas  of  the 
rolling  distant  countryside.  Such  a  landscape 
could  only  be  seen  by  ascending  local  promon- 
tories, and  few  locations  provide  this  opportun- 
ity. 

Stoneman  Lake,  a  caldera,  and  Apache  Maid 
Mountain  provide  outstanding  feature  land- 
scapes in  the  study  area.  Distinct,  visually  im- 
portant landscapes  are  relatively  few  on  this 
uniform  plateau.  , 

Significant  "macro"  landscapes  are  provided 
at  locations  where  one  can  look  out  across  the 
Verde  Valley  —  which  suggests  the  oppor- 
tunities for  enhancing  these  views  by  carefully 
planned  forest  clearing  along  certain  roadways. 

The  glades  surrounded  by  forest  referred  to 
by  Lt.  Beale,  may  be  what  are  now  known  locally 
as  "cienegas"  or  "parks."  Litton  would  call 
these  enclosed  landscapes:  they  are  defined  by 
their  "wall"  and  "floor"  characteristics.  Land 
managers  can  create  enclosed  landscapes  by 
cutting  harmonious  forest  openings.  Focal  land- 
scapes, on  the  other  hand,  can  occur  anywhere 
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that  roads,  landforms,  vegetative  patterns,  or 
waterways  lead  the  observer's  eye  to  a  point  of 
convergence. 

Canopied  landscapes,  characterized  by  im- 
mediate overhead  trees,  abound  for  the  hiker  in 
the  ponderosa  pine  forest  but  are  observed  by 
the  motorist  only  when  he  is  travelling  narrow, 
unimproved  roads. 

Detail  landscapes,  tied  to  minute  segments  of 
the  immediate  foreground,  are  an  important 
consideration  in  making  esthetic  evaluations  of 
forest  treatments.  It  is  the  "micro"  component 
of  the  landscape  which  includes  such  things  as 
the  bole  of  a  mature  ponderosa  pine  tree, 
lichen-covered  rocks,  or  bracken  fern  on  the 
forest  floor.  It  can  also  include  the  slash, 
stumps,  and  soil  disturbance  that  typify  some 
cutover  areas. 

Ephemeral  landscapes  exist  only  for  a  brief 
period.  Their  creation  and  disappearance  are 
generally  caused  by  climatic  or  phenological 
events  such  as  sunsets  and  autumnal  coloration. 
The  Gambel  oaks,  aspen,  and  occasional  maples 
associated  with  the  pine  forests  provide  an 
abundance  of  fall  colors  which  are  widely  ap- 
preciated in  Arizona.  Likewise  the  effect  of  man 
may  be  ephemeral,  depending  on  the  perma- 
nence of  his  impact. 

Other  concepts  of  landscapes  quality  dealing 
with  "variety,"  "deviations,"  and  "contrast" 
(U.S.  Forest  Service  1973)  also  have  implica- 
tions for  forest  treatments.  Since  they  imply 
changes  in  forest  cover,  however,  we  will  deal 
with  them  under  treatment  responses. 

Many  people  feel  that  quality  landscapes  are 
characterized  by  a  pleasing  degree  of  variety, 
with  the  various  components  of  the  landscape 
appearing  to  be  in  harmony  with  each  other. 
Harmony  demands  that  variety  be  of  a  compati- 
ble size  scale;  also,  it  is  important  that  small- 
scale  variations  should  not  be  repeated  indefi- 
nitely, nor  should  large-scale  variations  be 
easily  recognized.  These  concepts  are  obviously 
quite  subjective. 

The  overall  landscape  quality  of  this  pine 
area,  although  pleasing,  can  be  improved.  It  can 
also  be  damaged  by  inappropriate  land  man- 
agement practices.  Like  other  resources,  the 
scenic  resource  requires  sound  investments 
guided  by  skillful  management  in  order  to  attain 
its  potential  value,  especially  in  view  of  the  in- 
creasing and  competitive  demands  being  made 
on  the  forests  of  the  Rim  Country. 

Treatment  Tests 

Seven  treatments  are  presently  being 
evaluated  in  the  pine  type  on  Beaver  Creek.  Two 


more  are  scheduled  in  1974.  With  one  exception, 
a  wildfire,  all  treatments  are  on  calibrated 
watersheds.  All  are  timber  overstory  removal 
treatments  except  for  grazing. 

The  earlier  treatments  were  designed  to  in- 
crease streamflow,  but  were  evaluated  in  mul- 
tiple use  terms.  The  two  most  recent  treat- 
ments, on  Watersheds  14  and  16,  incorporated 
numerous  features  to  enhance  timber,  wildlife, 
and  esthetic  values.  The  two  treatments 
scheduled  in  1974  are  designed  primarily  to  en- 
hance wildlife  (Watershed  10)  and  timber 
(Watershed  8).  The  early  treatments  have  been 
described  by  Brown  (1970, 1971)  and  O'Connell 
and  Brown  (1972). 

Two  basic  kinds  of  timber  overstory  removal 
are  considered:  thinning  and  cleared  openmgs. 
Different  basal  areas  and  spacing  of  openings 
are  being  tested  to  evaluate  a  range  of  treat- 
ment intensities.  Clearing  all  the  woody  vegeta- 
tion is  being  tested  as  the  most  severe  level  of 
treatment  to  determine  maximum  possible 
water  yields. 

With  the  recent  shift  of  objectives  toward  de- 
velopment of  multiple  use  management  models, 
the  treatments  are  now  viewed  as  providing  val- 
idation data  for  mathematical  models  that  can 
be  used  to  estimate  treatment  response  for  a 
wide  range  of  conditions.  These  models  will  be 
discussed  in  a  later  section. 


Description  of  Treatments 

Thinning  by  group  selection  on  Watershed  17 
in  1969  (fig.  8). — Seventy-five  percent  of  the  ini- 
tial 120  ft"  basal  area  was  removed  by  thinning, 
leaving  evenaged  groups  with  average  basal 
areas  of  30  ft^/acre  (basal  areas  are  total,  for  all 
species).  All  Gambel  oaks  over  15  inches  d.b.h. 
except  den  trees  were  removed,  leaving  5  ft^/ 
acre  of  oak  basal  area;  all  alligator  junipers 
were  removed.  The  slash  was  windrowed. 
Theoretically,  the  treatment  would  be  repeated 
every  20  years. 

Stripcut  for  water  yield  on  Watershed  9  in 
1967-68  (fig.  9).— Thirty-two  percent  of  the 
watershed  was  clearcut  in  uniform  strips  60 
feet  wide,  and  slash  was  pushed  to  the  center  of 
cut  strips  and  partially  burned.  Spacing  be- 
tween cut  strips  is  120  feet.  Pretreatment  basal 
area,  all  species,  was  119  ft^/acre;  posttreat- 
ment  was  81  ft'.  Cut  strips  will  be  allowed  to 
restock  naturally.  The  area  between  cut  strips 
was  not  treated.  The  treatment  is  theoretically 
planned  to  be  repeated  every  40  years,  so  that  in 
a  rotation  of  120  years  the  whole  watershed 
would  be  cut. 
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Strip  shelterwood  cut  on  Watershed  14  in 

1970-  71  (fig.  10).— One-third  of  the  watershed 
was  cleared  in  irregular  cut  strips  averaging  60 
feet  wide.  Intervening  leave  strips  which  aver- 
age 120  feet  wide  were  thinned  to  80  ft-  basal 
area  by  a  silvicultural  cut  favoring  tree  diam- 
eter classes  that  are  in  short  supply  on  the 
Coconino  National  Forest.  The  original  basal 
area  of  all  species  was  121  ft-/acre;  the  post- 
treatment  basal  area  was  60  ft-  overall,  yielding 
an  overall  reduction  of  50  percent.  Heavy  con- 
centrations of  slash  from  leave  strips  were 
pushed  into  piles  in  the  cut  strips  and  burned. 
Pine  seedlings  were  planted  in  the  cut  strips  on 
better  sites.  Esthetic  areas  were  preserved 
from  any  treatment,  and  Gambel  oaks  less  than 
15  inches  d.b.h.  were  left  in  both  cut  and  leave 
strips  to  benefit  wildlife  and  esthetics.  Trees 
will  be  thinned  back  to  80  ft^/acre  at  20-year 
intervals,  and  nev..  strips  will  be  cut  at  40-year 
intervals. 

Strip  shelterwood  cut  on  Watershed  16  in 

1971-  72. — Fifty  percent  of  the  watershed  was 
cut  in  irregular  strips  averaging  60  feet  wide.  In 
the  intervening  leave  strips,  which  also  average 
60  feet  wide,  the  trees  were  thinned  to  80  ft^ 
basal  area.  Sixty-five  percent  of  the  basal  area 
was  removed  overall.  Gambel  oak  and  slash 
were  treated  as  in  Watershed  14.  Esthetic  areas 
were  preserved  from  any  treatment.  At  20-year 


Figure  9. — One-third  of  Watershed  9  was  cleared  in  60-foot  wide  strips,  with  no  treatment  of  the  inter- 
vening areas.  The  strip  width  was  chosen  to  maximize  trapping  and  retention  of  snow.  The  strips  were 
oriented  to  lead  water  directly  into  the  streams. 


Figure  8. — In  1 969,  75  percent  of  the  basal  area  was 
removed  on  Watershed  17  in  a  heavy  thinning 
treatment. 
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intervals  the  stands  will  be  thinned  back  to  80 
ft^/acre  of  basal  area,  and  at  60-year  intervals 
the  alternate  strips  will  be  cut. 

Total  clearcut  on  Watershed  12  in  1966-67  (fig. 

11)  . — This  watershed  was  essentially  removed 
from  timber  production.  The  slash  was  wind- 
rowed  and  oak  sprouts  were  partially  controlled 
with  chemical  sprays.  While  some  natural 
restocking  is  taking  place,  it  will  be  many  years 
before  any  timber-oriented  work  will  be  done 
unless  the  area  is  planted.  This  treatment  is  not 
intended  to  have  operational  potential  but 
rather  is  an  analytical  benchmark  against  which 
to  compare  less  severe  treatments. 

Grazing  on  converted  Watershed  11  (fig. 

12)  . — Sixty  percent  of  the  perennial  grass  on  the 
watershed  has  been  utilized  by  cattle  grazing  in 
the  spring  and  fall,  starting  in  the  fall  of  1967. 
The  watershed  had  been  converted  from  pine  to 
grass  in  1958. 

Wildfire  of  severe  and  moderate  intensities  in 
1972. — The  effects  of  wildfire  are  being  studied 
on  two  small  uncalibrated  watersheds,  one 
burned  severely  and  the  other  moderately.  The 
Rattle  Burn  fire  occurred  in  the  spring  of  1972, 
30  miles  southwest  of  Flagstaff.  Shortly  after- 
wards three  small  watersheds  were  in- 
strumented, two  on  the  burn  and  one  nearby  for 
a  control.  Special  emphasis  in  the  evaluation  of 
this  treatment  is  being  placed  on  wildlife  ef- 
fects, nutrient  movement,  and  water  infiltration 
into  the  soil. 


Treatment  Responses 

Streamflow  responses  have  been  evaluated 
on  six  Watersheds:  12  —  clearcut;  9  —  one-third 
regular  stripcut;  17 —  three-fourths  thin;  14 — 
one-third  irregular  stripcut  and  thin;  16 —  one- 
half  irregular  stripcut  and  thin;  and  11 —  graz- 
ing. Total  sediment  yield  has  been  evaluated 
only  on  Watersheds  12,  9,  17,  14,  and  16.  An 
unusually  large  storm  in  September  1970  pro- 
vided an  additional  basis  for  evaluating  flood 
response  on  Watersheds  12,  8, 17,  and  11.  Inten- 
sive suspended  sediment  sampling  provides  the 
only  basis  for  evaluating  sediment  from  the 
grazing  treatment  on  Watershed  11.  Herbage 
responses  have  been  intensively  evaluated  on 
Watersheds  12,  9,  and  11.  Herbage  production 
models,  now  under  development,  will  be  used 
later  to  predict  responses  to  the  other  treat- 
ments. 

Streamflow. — Responses  of  winter  stream- 
flow  to  the  vegetative  treatments  on  the  pine 
watersheds  are  summarized  in  table  2.  Mean 
winter  streamflow  for  the  untreated  condition 
is  based  on  variable  periods  of  record,  and  has 
been  adjusted  by  covariance  to  minimize  uncon- 
trolled variation  due  to  climate  and  other  fac- 
tors not  associated  with  treatment  (Kovner  1956, 
Kovner  and  Evans  1954,  Wilm  1943).  This 
analysis  is  designed  to  detect  treatment  re- 
sponses of  greater  than  15  percent  at  appropri- 
ate risk  levels.  Responses  for  individual  years 
following  treatment  are  determined  by  sub- 
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Figure  10. — One-third  of  Watershed  14  was  cleared  in  1971  in  a  pattern  of  irregular  strips  which  vary  in 
width  from  30  to  90  feet.  The  intervening  strips  have  been  thinned  to  80  ft^  basal  area  to  improve  timber 
growth.  The  cut  strips  on  better  sites  have  been  planted  to  assure  continued  timber  production.  Oak 
trees  are  left  standing  in  the  cut  strips  to  enhance  esthetic  and  wildlife  values.  Three  especially  attrac- 
tive "esthetic  areas"  were  not  treated. 
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Figure  11. —  Aerial  view  of  pine  clearcut  on 
Watershed  12.  The  slash  was  windrowed  in  a  pat- 
tern to  trap  snow  and  expedite  water  movement 
into  the  streams. 


Table  2. --Summary  of  effects  of  treatments  on  winter  streamflow  from  the  pine  watersheds 


1,^      ...^      ^      ^           ,,  ^      ^    ,  Difference  between  actual  and        ^           _          4.  i       i  „f 

Watershed,  by  treatment,          Untreated  ,.       ,                     ,                    Mean      Percent  Level  or 

.   ^  predicted  streamflow  by  years       .-rr        a-cc  ^      •  -f 

and  vear                        mean  winter  *^      r  ^^     ■      ^      ^      I2              differ-  differ-  signifi- 

^        r,     1  following  treatment 

treatment  completed  streamflow   ^3  3d   ITEIt  — 5th  &tF  ^""^^  cance 


Watershed  12-- 

100%  clearcut  1967 

Watershed 

32%  clearcut  in 

uniform  strips  I968 

Watershed  11-- 

()0%  ut  i  1  i  zat  ion  of 
perennial  grass  1967 

Watershed  I7-- 

75%  thin  1969 

Watershed  l^*-- 

33%  pine  clearcut  in 
irregular  strips, 
thinning  between 
strips;  50%  basal 
area  removed  overall  1970 

Watershed  I6-- 
50%  clearcut  in 
irregular  strips, 
thinning  between 
strips;  65%  basal 
area  removed  overall  1972 


6.04  3.79  0.92  1.81  1.47 

6.70  1.98  .61  .34  .8't 

3.16  .3^  .24  .11  .29 

7.63  .85  1.45  1.51  2.93 

4.71  .71  .70  1.61 
5.45  5.60 


1.39    3.29    2.03  34  0.01 

1.74             1.10  16  .05 

-.20      .78      .26  8  (^) 

1.68  22  .05 

1.01  21  .05 

5.60  103  .01 


^  Adjusted  by  covariance  for  years  1960-73  for  watersheds  12,  9,   14;   1961-73,  watershed  11;  1963-73, 
watersheds  16  and  17. 

^  Predicted  streamflow  is  that  which  would  be  expected  to  occur   in  a  given  year  had  the  watershed 

not  been  treated. 
^  Not  significant  at  0.05  level. 
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trading  the  flow  that  would  be  predicted  to 
occur,  assuming  pretreatment  vegetative  con- 
ditions, from  the  actual.  Average  treatment  re- 
sponses are  expressed  quantitatively  in  inches 
and  as  a  percent  of  the  pretreatment  mean. 

The  largest  response  in  winter  streamflow 
was  from  an  irregular  stripcut  and  thinning 
treatment  on  Watershed  16  which  removed  65 
percent  of  the  timber  basal  area.  The  response 
averaged  5.6  inches  or  103  percent  for  the  first 
year  after  treatment,  and  was  significant  at  the 
1  percent  level.  This  very  large  response  should 
be  interpreted  with  caution.  It  is  partially  the 
combined  result  of  the  unusually  wet  1973 
winter  and  the  large  increase  that  is  usual  the 
first  posttreatment  year.  The  theory  behind  the 
irregular  stripcut  and  thin  treatment  is  that  de- 
creased water  loss  by  transpiration  and  in- 
creased efficiency  of  runoff  transport  in  the  cut 
strips  will  increase  water  yield  from  the 
watershed.  More  efficient  waterways  for 
transporting  surface  runoff  into  the  channels 
are  believed  to  be  formed  by  orienting  cut  strips 
perpendicular  to  the  slopes  and  by  reducing  the 
soil  water  deficit  through  tree  removal.  More 
snow  should  also  reach  the  ground  in  the  cut 
strips  due  to  reductions  in  interception  losses 
and  in  redistribution  of  snow  by  the  wind.  The 
snowpack  in  the  east-west  oriented  strips  is  par- 
tially shaded,  thus  reducing  evaporation  losses 
during  the  winter.  This  also  produces  more  effi- 
cient runoff  by  delaying  melt  until  the  ambient 
temperature  is  relatively  high. 

The  second  largest  response  in  winter 
streamflow  came  from  the  100  percent  clearcut 
on  Watershed  12.  The  runoff  increase  averaged 
2.0  inches  per  year  or  34  percent  over  a  6-year 
period,  and  was  significant  at  the  1  percent 
level.  The  increased  water  yield  from  this 
clearcut  treatment  is  primarily  due  to  a  de- 
crease in  water  loss  by  transpiration.  By  reduc- 
ing the  soil  moisture  deficit  on  the  watershed, 
more  runoff  is  realized  from  the  melting  snow- 
pack.  In  addition,  the  windrowed  slash  traps 
snow  and  delays  melt  rates  on  the  lee  side.  With 
delayed  snowmelt,  snow  remains  on  the 
watershed  until  the  ambient  temperature  rises 
high  enough  to  melt  the  snow  rapidly,  and  more 
surface  runoff  thereby  reaches  the  channels. 
The  efficiency  of  surface  runoff  transport  is 
also  increased  because  the  windrows,  which  are 
directed  downslope,  shade  the  soil  on  their  lee 
side.  This  reduces  soil  moisture  loss  and  pro- 
vides more  efficient  waterways  for  transport- 
ing surface  runoff  into  the  channels. 

The  third  largest  response  resulted  from  the 
removal  of  75  percent  of  the  timber  basal  area 
by  thinning  on  Watershed  17.  The  average  in- 
crease was  1.7  inches  per  year  or  22  percent  for 
a  4-year  period,  and  was  significant  at  the  5  per- 


cent level.  Again  the  response  is  apparently  due 
largely  to  greater  efficiency  in  surface  runoff 
and  reduced  transpiration  loss.  The  residual 
slash  was  piled  into  windrows  and  had  an  influ- 
ence on  the  snowpack  similar  to  that  on  the 
clearcut  Watershed  12. 

Watershed  14,  which  was  stripcut  and 
thinned,  provided  a  1.0  inch  per  year  or  21  per- 
cent response  over  a  3-year  period,  significant 
at  the  5  percent  level.  Watershed  14  can  be  ex- 
pected to  behave  similarly  to  Watershed  16  but 
probably  not  as  efficiently  because  only  a  total 
of  50  percent  of  the  timber  basal  area  was 
removed. 

Watershed  9  had  32  percent  of  its  timber  re- 
moved in  uniform  downslope  strips.  Stream- 
flow  increase  was  1.1  inches  per  year  or  16  per- 
cent over  a  5-year  period,  significant  at  the  5 
percent  level.  Watershed  9  is  expected  to  be- 
have similarly  to  Watersheds  14  and  16.  The 
vegetation  treatment  differs  on  Watershed  9  in 
that  the  cut  strips  are  uniform  in  width,  the 
leave  strips  were  not  thinned,  and  no  oak  trees 
were  left  in  the  cut  strips. 

Watershed  11,  which  has  been  grazed  after 
conversion  to  a  grass  cover,  shows  only  an  8 
percent  increase  after  a  6-year  period.  This  in- 
crease is  not  statistically  significant  at  the  5 
percent  level,  and  does  not  approach  practical 
importance.  It  was  theorized  when  this  treat- 
ment was  planned  that  grazing  might  reduce 
herbaceous  plant  growth  and  result  in  an  in- 
crease in  water  yield  due  to  reduced  transpira- 
tion losses.  It  was  also  thought  that  compaction 
of  the  soil  by  the  cattle  might  contribute  to  the 
efficiency  of  surface  runoff.  Plant  growth  has 
changed  negligibly,  however,  (see  table  5),  and 
soil  compaction  apparently  was  not  sufficient  to 
produce  an  important  increase  in  surface 
runoff. 

A  special  commentary  is  in  order  regarding 
the  unusually  wet  winter  of  1972-73,  and  its  ef- 
fect on  watershed  treatment  responses.  The 
highest  previously  recorded  winter  precipita- 
tion for  pine  watersheds  on  Beaver  Creek  was  in 
water  year  1965,  and  ranged  from  21  to  28  in- 
ches. The  resulting  runoff  was  also  the  previous 
high  on  record,  ranging  from  6.3  to  14.1  inches. 
In  the  winter  of  1972-73  these  watersheds  re- 
ceived between  32  and  41  inches  or  an  average 
of  1.5  times  more  precipitation  than  in  1965. 
Runoff  was  11.1  to  27.8  inches  or  an  average  of 
1.8  times  more  than  the  previous  high.  Re- 
sponses to  treatment  during  water  year  1973 
exceeded  all  but  the  first-year  responses  on 
Watersheds  9  and  12  (table  2). 

Flood  peak. — The  major  storm  which  struck 
Beaver  Creek  in  September  1970  provided  an 
opportunity  to  assess  the  effect  of  the  pine 
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treatments  on  flood  peaks."  Estimates  of  treat- 
ment effects  are  summarized  in  table  3.  These 
responses  were  determined  graphically,  based 
on  the  relationship  between  peak  discharge  and 
60-minute  precipitation  intensity  (fig.  13).  The 
assumption  here  is  that  the  departures  of  a 
treated  watershed  from  the  untreated 
watershed  relationship  represent  a  treatment 
response.  Table  3  also  contains  the  estimated 
recurrence  intervals  for  the  observed  peak  dis- 
charges. 

Peak  discharge  increased  most  —  750  f t^/s/mi^ 
—  on  Watershed  12.  The  second  largest  re- 
sponse was  400  ftVs/mi^  from  Watershed  11, 
which  had  been  clearcut  and  seeded  to  grass  in 
1958,  12  years  prior  to  the  storm. 

Responses  which  were  apparently  due  to  the 
complete  or  large  reduction  in  timber  and  as- 
sociated litter  cover  from  Watersheds  11,  12, 
and  17  ranged  from  350  to  750  ft^/s/mi%  or  about 
twice  the  peak  discharge  expected  from  an  un- 

"Bafcer,  Malchus  B.,Jr.,  Harry  E.  Brown,  and  Norman  E. 
Champagne,  Jr.  Hydrologic  performance  of  the  Beaver 
Creek  watersheds  during  a  100-year  storm.  Paper  pre- 
sented at  Amer.  Geophys.  Union  meeting,  San  Francisco, 
Calif,  19  p.  Dec.  1971. 


treated  watershed  with  similar  precipitation  in- 
tensity. On  the  average,  92  percent  of  the  timber 
basal  area  had  been  removed  from  these  water- 
sheds. They  were  subjected  to  5.6  inches  total 
rainfall  with  1.1  inches  maximum  60-minute  in- 
tensity. 

The  recurrence  interval  for  the  flood  runoff 
ranged  between  100  and  200  years  for  the 
treated  pine  watersheds.  This  corresponds 
favorably  with  the  recurrence  interval  of  100 
years  for  24-hour  rainfall  of  5.5  to  6.0  inches 
(U.S.  Department  of  Commerce  1967). 

Sediment. — Annual  sediment  yields  from 
clearcut  Watershed  12  have  varied  from  0.01  to 
27  tons  per  acre  since  treatment.  On  stripcut 
Watershed  9,  with  a  residual  basal  area  of  81 
ft'^/acre,  they  have  varied  from  0.01  to  1.3  tons 
per  acre.  On  heavily  thinned  Watershed  17,  with 
30  ft^  basal  area,  they  have  been  from  0.03  to 
0.32  ton  per  acre. 

Watershed  11  does  not  have  an  installation  to 
measure  total  sediment  yield,  but  suspended 
sediment  samples  were  collected  periodically 
before  and  after  the  grazing  treatment  was  im- 
posed. Average  suspended  sediment  concentra- 


Table  3. --Summary  of  effects  of  treatments  on  flood  peaks  for  the  September  1970  storm  on  the  pine 
watersheds   (all  discharges  rounded  to  nearest  50  ft^/s/mi^) 


Flood  peak 


Watershed,  by  treatment, 
and  year 
treatment  completed 


Est  i  mated , 
unt  reated 
condit  ions 


Actual 


Est  i  mated 
d  i  f f e rence 
due  to 
t  reatment 


Estimated  Since 
repurrence  treat- 
i  nterval ^  ment 


Watershed  12-- 
100^  clearcut 


  ft^/s/mi^ 

1967  hSO  1,200 


750 


-  Years 
150-200 


Watershed  9"- 
32^  clearcut  in 

uniform  strips  1968  500  600  '^\00  100-150  2 

Watershed  11-- 

60^  ut  i 1 i  zat  ion  of 

perennial  grass  1967  ^50  850  hOO  100-150  ^12,(3) 

Watershed  17"- 

75^  thin  1969  'too  750  350  100-150  1 

Watershed  l^t-- 

33%  pine  clearcut  in 
irregular  strips, 
thinning  between 
strips;  50%  basal 

area  removed  overall   I970  450  550  ^joo  IOO-I5O  C*) 


Estimated  recurrence  intervals  were  determined  by  the  U.S.  Geological  Survey,  Tucson,  Arizona, 
using  regionalized  peak  discharge  data  from  the  Beaver  Creek  watersheds. 
Only  a  marginal  response  if  indeed  significant. 

Twelve  years  since  clearcut  and  3  years  since  grazing  treatment  imposed. 
Treatment  'tO  percent  completed. 
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Figure  13. — 
Flood  peak  response 
determined  by  relationship 
of  peak  discharge  to  maximum 
60-minute  rainfall  intensity. 


tion  prior  to  grazing  was  176  p/m  (38  samples) 
and  the  average  following  treatment  was  113 
p/m  (29  samples).  The  difference  between  these 
two  averages  is  not  significant  at  the  20  percent 
level. 

The  largest  sediment  yield,  27  tons  per  acre, 
was  produced  on  Watershed  12  after  it  had  just 
been  clearcut.  Because  all  slash  had  been 
pushed  into  windrows,  nearly  100  percent  of  the 
surface  soil  was  disturbed  by  the  piling  equip- 
ment. The  entire  treatment  was  completed  in 
June  1967.  On  July  31,  1967,  this  watershed  re- 
ceived 2.2  inches  of  precipitation  with  a  max- 
imum 30-minute  intensity  of  2.8  inches  per  hour. 
The  resulting  2.1  inches  of  runoff  overtopped 
the  flume  and  washed  out  the  sediment  basin 
below  it.  The  sediment  yield  of  27  tons  per  acre 
was  estimated  from  measurements  in  the  field 
and  substantiated  by  measurements  obtained 
from  pre-and  poststorm  aerial  photographs  of 
the  watershed. 

From  the  sediment  yield  records  obtained 
during  the  past  11  years  on  Beaver  Creek,  and 
from  knowledge  of  the  sediment  losses  result- 
ing from  various  treatment  intensities  and 
storm  frequencies,  it  appears  that  27  tons  per 
acre  is  approaching  the  maximum  sediment  loss 
potential  from  watersheds  with  similar  charac- 
teristics in  the  ponderosa  pine  type. 

It  is  difficult  to  develop  a  relationship  to  ade- 
quately predict  sediment  yields  because  of  the 
limited  amount  of  data.  However,  an  initial  at- 
tempt will  be  described  in  a  later  section,  using 
maximum  annual  30-minute  precipitation  in- 
tensity and  timber  basal  area  as  independent 
variables. 


Water  quality.— TDS  (total  dissolved  solids), 
which  are  indicative  of  overall  water  quality, 
were  apparently  not  affected  by  treatments.  Al- 
though the  mean  TDS  from  treated  watersheds 
was  nearly  5  milligrams  per  liter  greater  than 
from  untreated  watersheds,  a  difference  of 
nearly  10  percent,  analyses  were  not  statisti- 
cally conclusive. 

Analysis  of  correlation  between  discharge  in 
the  flume  at  the  time  of  sample  collection  and 
TDS  in  the  samples  from  the  pine  watersheds 
yielded  a  coefficient  of  -0.44,  not  significant  at 
the  5  percent  level.  Thus  no  clear  relation  is 
evident  between  rate  of  flow  and  quality  of 
runoff  water. 

Samples  from  both  treated  and  untreated 
watersheds  met  water  quality  standards  for 
drinking  water,  fish  and  aquatic  life,  and  unre- 
stricted use  for  irrigation.  Only  iron  ap- 
proached an  unacceptable  level  for  drinking 
water.  All  samples  had  extremely  low  sodium 
absorption  ratios.  Electrical  conductivity  val- 
ues of  all  pine  watershed  samples  were  72  /x 
mho/cm  or  under. 

Based  on  the  classification  of  Durfor  and 
Becker  (1964,  p.  27),  water  samples  from  the 
pine  watersheds  were  soft  (0  to  60  mg/1  of 
CaC03). 

Suspended  sediment  was  determined  on  sam- 
ples taken  concurrently  with  those  reported 
here.  Correlation  analysis  showed  no  consistent 
relation  between  suspended  sediment  and  TDS 
of  the  waters. 

Timber. — No  direct  measures  of  timber  re- 
sponses have  been  possible  because  the  water- 
sheds have  not  had  time  to  respond  to  the  treat- 
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ments.  Thus  the  only  means  available  for  de- 
termining response  has  been  the  simulation 
model  described  in  a  later  section.  This  model 
uses  Beaver  Creek  inventory  data  consisting  of 
an  average  stand  table  for  each  soil  type.  Esti- 
mates for  rate  of  natural  regeneration  and  mor- 
tality are  determined  for  each  soil  type.  Func- 
tions are  computed  for  growth  and  the  effect  of 
density  on  growth. 

All  of  the  treatments  were  simulated  for  a  full 
rotation,  and  the  total  products  harvested  are 
estimated  in  table  4.  Past  management  was  also 
simulated  for  each  watershed,  and  these  results 
are  also  included  in  the  table.  Past  management 
consisted  of  harvest  cuts  every  40  years  to  re- 
move overmature  and  dying  trees,  salvageable 
snags,  and  some  mature  and  older  blackjacks. 
No  timber  stand  improvement  work  was  done, 
and  no  pulpwood  was  sold. 

Volumes  removed  in  the  initial  treatment  cuts 
were  omitted  from  the  simulation  analysis  in 
order  to  best  describe  the  long-term  effects  of 
the  treatments  without  the  bias  introduced  by 
the  initial  cut.  Consequently  the  cleared  Water- 


sheds, 11  and  12,  produce  no  timber  under  the 
alternative  treatment  prescriptions,  and  a  net 
loss  in  timber  products  is  shown.  Every  40 
years,  one-third  of  Watershed  9  would  be 
cleared,  resulting  in  young  stands  capable  of 
rapid  cubic-foot  growth  until  tree  growth  di- 
minishes due  to  increasing  stand  density.  This 
pattern  causes  rapid  growth  on  young  trees  but 
slow  growth  in  older  saw  log-sized  stands,  with  a 
net  gain  in  cubic-foot  volume  but  net  loss  in 
board-foot  volume.  While  one-third  of 
Watershed  14  is  cut  every  40  years,  timber  pro- 
duction is  higher  than  under  past  management 
for  two  reasons:  the  cut  strips  are  planted,  thus 
keeping  the  watershed  under  full  production, 
and  the  leave  strips  are  thinned  for  stand  im- 
provement and  growth.  Estimates  for 
Watershed  17  indicate  that  it  may  be  possible  to 
maintain  good  volume  production  at  very  low 
stocking  levels  by  intensive  management. 

Range. — With  the  exception  of  the  grazing 
treatment  on  Watershed  11,  the  Beaver  Creek 
treatments  consist  of  reductions  of  timber 


Table  4. --Pred i cted  treatment  effects  on  wood  products  sold  from  the  pine  watersheds  for  a  120-yr 

rotat  i  on 


Watershed,  by  treatment, 

Volume  per 

acre 

and  year 
treatment  completed 

Under  past 

Under  alternative-treatment 
management^  prescriptions^ 

D  i  f f erence 

ft' 

fbm 

ft' 

fbm 

ft' 

fbm 

Watershed  12-- 

100^  clearcut  1 

967 

800 

4,200 

0 

0 

-800 

-4,200 

Watershed  9-- 

32Z  clearcut  in 

un  i  form  st  r  i  ps  1 

968 

1  ,500 

7,500 

2,000 

6,300 

500 

-1 ,200 

Watershed  11-- 

60%  ut  i 1 izat  ion  of 

perennial  grass  1 

967 

500 

2,800 

0 

0 

-500 

-2,800 

Watershed  1?-- 

7S%  thin  1 

969 

1  ,400 

6,800 

3,900 

1 1  ,700 

2,500 

4,900 

Watershed  14-- 

33%  pine  clearcut  in 

irregular  strips, 

thinning  between 

strips;  50%  basal 

area  removed  overal 1 

970 

1  ,100 

5,500 

3,300 

1 1  ,200 

2,200 

5,700 

^Past  management  consisted  of 

harvest  cuts  every  40  years 

to  remove 

overmature 

and 

dy  i  ng 

trees , 

salvageable  snags,  and 

some 

mature  and  o 

1 der  bl ackj acks . 

^Volumes  removed   in  the 

i  n  i  t  i 

al  cut  were 

omitted   in  order 

to  best  describe  the 

1  ong 

-term 

effects 

of  the  treatment  without  the 

bias  introduced  by  the  initial  cut. 
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overstory  by  differing  degrees  and  in  different 
patterns.  Such  reductions  invariably  result  in 
an  increase  in  the  production  of  herbaceous 
plants,  and  often  of  shrubby  understory  plants 
as  well.  These  responses  increase  with  increas- 
ing degrees  of  overstory  removal. 

On  Watershed  12,  all  timber  overstory  was 
removed.  The  native  herbaceous  and  shrubby 
plants  increased  an  average  of  504  pounds  per 
acre  for  the  five  posttreatment  years  (table  5). 
This  increase  in  understory  plant  yields  means 
greater  livestock  and  big-game  carrying  capa- 
cities. If  the  cattle  carrying  capacities  are  based 
on:  (1)  grass  production  only,  (2)  assumption 
that  750  pounds  of  allowable  forage  disappear- 
ance represents  one  animal-unit  of  carrying 
capacity,  and  (3)  40  percent  maximum  allow- 


able use,  then  the  carrying  capacity  for  cattle 
would  have  increased  from  0.11  animal-unit 
month  per  acre  to  0.23  AUM  per  acre  as  a  result 
of  treatment. 

Watershed  9  was  stripcut  with  no  thinning  of 
the  leave  strips,  thus  herbage  response  can  only 
be  expected  on  one-third  of  the  area.  By  the 
second  year  after  treatment  the  average  in- 
crease in  production  of  understory  plants  was 
statistically  significant,  but  not  yet  of  sufficient 
magnitude  to  have  practical  importance.  The 
sampling  of  the  stripcut  treatment  was  rela- 
tively intense,  but  limited  to  specific  combina- 
tions of  slope,  aspect,  and  soils.  Thus  the  results 
do  not  represent  a  watershed  mean. 

The  treatment  of  Watershed  11  consisted  of 
spring  and  fall  grazing  by  livestock,  with  the 


Table  5. --Summary  of  effects  of  treatments  on  herbaceous  and  shrubby  plant  production  for  the  pine 

watersheds  on  Beaver  Creek 


_     _        Difference  between  actual  and  -  ... 

Watershed,  by  treatment,  ,  predicted  production  by  years  statistical 

and  year  Plant  group      Auction  following  treatment^  significance 

treatment  completed  '      untreated   ^ver-      °^  ^^'^f 

conditions     1st      2d      3d  'tth  5th  change"^ 
 age  ^ 

-  -  -  -  Pounds  per  acre  -  -  -  - 


Watershed  12-- 

Grasses 

200 

176 

252 

17'* 

262 

319 

237 

0.01 

100^  clearcut 

1967 

Forbs , ha  1 f- 

shrubs 

109 

313 

168 

156 

21  1 

382 

2k6 

.02 

Sh  rubs 

20 

-20 

1 

1  1 

66 

21 

NS 

Watershed  9"-- 

Grasses 

30 

12 

10 

1  1 

.05 

32%  clearcut  in 

Forbs , ha  1 f- 

shrubs 

26 

32 

.01 

uniform  strips 

1968 

Shrubs 

1  1 

3 

12 

8 

NS 

Watershed  11-- 

Grasses 

390 

-118 

-1 

215 

303 

-12 

77 

NS 

GOZ  ut  i 1 i  zat  ion  of 

Forbs , ha  1 f- 

sh  rubs 

339 

-131 

-104 

-181 

-61 

-5 

-96 

.10 

perennial  grass 

1967 

Shrubs 

0 

1 

2 

0 

1 

0 

1 

NS 

Watershed  17--                            Grasses  86 

75%  thin  1969      Forbs, half-shrubs  ^5 

Shrubs  15 

Watershed  l^t--                            Grasses  137 

33%  pine  clearcut  in               Forbs , ha  1 f-shrubs  86 

irregular  strips,                    Shrubs  2 
thinning  between 
strips;  50%  basal 

area  removed  overall  1970 


Predicted  production   is  an  estimate  of  production  that  would  be  expected  assuming  the  watershed 
had  not  been  treated. 

Different  sample  designs  and  sampling  intensities  were  used  on  the  different  watersheds;  therefore, 
the  reliabilities  of  the  estimates  differ. 

In  1968  and  1969,  attempts  were  made  to  control   Gambel  oak  sprouts  with  chemicals.     The  original 
treatment  plus  followup  chemical  application  suppressed  shrub   (browse)  production  for  several  years. 
Since  that  time,  Gambel  oak  sprout  production  has  been  steadily  increasing.     Additional  sampling 
has  indicated  that  oak  sprout  production  has  increased  three  to  five  times  since  treatment. 
Sampling  was   limited  to  sites  exceeding  10  percent  slope.     These  sites  have  greater  timber  densities 
and  therefore  low  initial   herbage  production. 
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intent  of  removing  one-half  of  the  standing  crop 
of  perennial  grasses  during  each  grazing 
period.  The  proportion  of  the  total  perennial 
grass  production  removed  yearly  was  60  per- 
cent for  the  first  4  years  of  treatment,  but 
dropped  considerably  the  fifth  year.  This  treat- 
ment caused  no  significant  change  in  the  pro- 
duction of  grasses,  but  did  apparently  reduce 
the  production  of  forbs.  If  we  again  base  the 
carrying  capacity  for  cattle  upon  the  grasses 
alone,  we  would  conclude  that  there  has  been  no 
change  in  carrying  capacity  as  a  result  of  this 
treatment. 

The  understory  plant  responses  on  Water- 
sheds 17  (3/4  thin)  and  14  (stripcut  plus  thin) 
have  not  been  directly  sampled.  It  appears  that 
after  several  posttreatment  years  these  water- 
sheds are  producing  about  100  pounds  of  addi- 
tional herbage  per  acre,  about  one-half  being 
grass.  The  dry  years  since  treatment  have  obvi- 
ously slowed  the  plant  responses.  Therefore, 
additional  recovery  time  is  needed  before  live- 
stock carrying  capacity  will  increase  meaning- 
fully. 

In  evaluating  the  effect  of  a  treatment  on 
livestock  carrying  capacity,  we  must  consider 
whether  the  animals  can  or  will  take  advantage 
of  the  increased  forage  yields.  In  the  case  of 
Watershed  12  the  topography  is  relatively  gen- 
tle, abundant  water  is  nearby,  and  the  cattle 
have  demonstrated  that  they  will  distribute 
themselves  over  the  watershed.  Here  we  can 
conclude  that  the  apparent  increase  in  carrying 
capacity  is  realistic.  On  other  areas  such  as 
Watershed  9,  any  increase  in  apparent  carrying 
capacity  will  have  to  be  tempered  by  the  knowl- 
edge that  livestock  will  probably  spend  little 
time  on  the  steeper  areas. 

Wildlife. — The  impact  of  watershed  treat- 
ments upon  wildlife  values  has  been  measured 
in  two  ways:  (1)  production  of  forage  plants 
preferred  by  deer,  and  (2)  big-game  fecal  drop- 
pings as  an  indirect  index  of  animal  use. 

We  do  not  know  what  a  given  increase  in  pre- 
ferred plants  means  in  terms  of  big-game  carry- 
ing capacity,  but  an  increase  in  the  forage  sup- 
ply should  be  beneficial  if  cover  requirements 
can  still  be  met  (Reynolds  1969).  Big-game  on 
Beaver  Creek  are  known  to  prefer  areas  with 
better-than-average  forage  supplies  (Clary  and 
Larson  1971).  The  overstory  reduction  treat- 
ments applied  to  Watersheds  12  and  9  resulted 
in  increased  production  of  the  preferred  deer 
forage  plants  (table  6).  The  livestock  grazing 
treatment  on  Watershed  11,  on  the  other  hand, 
reduced  preferred  plants. 

The  response  of  deer  and  elk  to  the  treat- 
ments, as  indexed  by  fecal  droppings,  roughly 
parallels  the  production  of  preferred  forage 


plants.  On  Watershed  12  with  the  highest  in- 
crease in  preferred  plants,  deer  use  increased 
noticeably  (table  7).  Big-game  use  changed  little 
on  Watershed  9,  which  showed  only  a  small  in- 
crease in  preferred  plants.  On  Watershed  11 
where  production  of  preferred  plants  has  been 
severely  depressed,  deer  use  remains  low, 
while  elk  use  has  plummeted.  On  the  other  hand 
both  deer  and  elk  appear  to  be  responding  to  the 
treatment  of  Watershed  17,  even  though  there 
does  not  appear  to  be  a  large  increase  in  herbage 
yields. 

In  general,  opening  the  ponderosa  pine 
canopy  to  enhance  the  understory  forage  supply 
improves  deer  and  elk  habitat  (Neff  1972).  It  is 
still  not  possible,  however,  to  predict  the  exact 
game  response  to  any  particular  treatment. 

The  question  of  quantifying  wildlife  cover  is  a 
difficult  one.  Most  wildlife  biologists  and  wild- 
life managers  avoid  attempting  to  quantify 
cover  requirements.  Clearcutting  Watershed  12 
(455  acres)  was  expected  by  some  to  create  a 
wasteland,  from  a  deer  standpoint.  In  fact,  how- 
ever, deer  use  of  this  watershed  is  several  times 
higher  than  before  treatment.  The  deer  are  ap- 
parently attracted  by  the  abundant  forage  sup- 
ply, while  the  slash  windrows  and  Gambel  oak 
sprouts  provide  sufficient  cover.  A  critical 
evaluation  of  the  effect  of  these  treatments  on 
wildlife  cover  requires  additional  understand- 
ing of  the  specific  requirements  of  individual 
species. 

Preliminary  predictive  models  based  on  pel- 
let counts  were  developed  from  several  differ- 
ent studies.  The  graphs  (fig.  14)  are  presented 
as  an  initial  basis  for  evaluating  the  impact  of  a 
given  treatment  on  deer  use.  The  three  differ- 
ent situations  illustrated  describe  the  relation 
of  deer  use  to  basal  area. 

The  shape  of  relationship  A  applies  to  cleared 
forest  openings  smaller  than  40  acres  where  the 
lack  of  cover  does  not  inhibit  deer  use  (Reynolds 
1969).  On  these  small  areas,  average  big-game 
use  will  likely  be  proportional  to  forage  supply 
(Reynolds  1962;  Patton  1969)  —  highest  at  low 
timber  basal  areas. 

Relationship  B  applies  to  areas  40  to  500  acres 
in  size,  surrounded  by  forest.  Under  these  cir- 
cumstances deer  use  peaked  on  a  heavily 
thinned  watershed  of  300  acres  where  there  was 
a  moderate  increase  in  forage  and  some  remain- 
ing tree  and  slash  cover  (Neff  1973).  The  shape 
of  relationship  C  is  applicable  to  broad,  un- 
iformly forested  areas  (Reynolds  1969)  where 
both  cover  and  forage  are  provided  on  the  same 
site. 

Future  intensive  land  management  will  un- 
doubtedly be  based  on  diversified  sites,  habitat 
types,  and  land  response  units,  resulting  in  a 
mosaic  of  treatment  patterns.  Under  these  cir- 
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Table  6. --Summary  of  effects  of  treatments  on  forage  plants  preferred  by 
deer  in  the  pine  type  on  Beaver  Creek 


Watershed,  by  treatment, 
and  year 
treatment  completed 


Deer  forage 
product  i  on , 
unt  reated 
cond  i  t  i  ons 


Difference  between  actual  and 
predicted  production  by  years 
 following  treatment   


1st 


2d 


3d     hth  5th 


Aver- 
age 


Stat  i  St  i  ca 1 
s  i  gn  i  f  i  cance 
of  average 
change 


Watershed  I2-- 

100%  clearcut  1967 

Watershed  3-- 
32%  clearcut  in 
uniform  strips  1968 

Watershed  1 1 -- 

60%  ut  i 1 i  zat  ion  of 
perennial  grass  1967 

Watershed  17"- 

75%  thin  1969 

Watershed  14-- 

33%  pine  clearcut  in 
irregular  strips, 
thinning  between 
strips;  50%  basal 
area  removed  overall  1970 


69 


16 


190 


30 


37 


Pounds  per  acre  - 
233     75  '*7 


21 


13 


86  166  121 


17 


•166    -63    -87    -80  -46 


0.05 


,01 


.10 


^Predicted  production  is  an  estimate  of  production  that  would  be  expected,  assuming  the  watershed 
had  not  been  treated. 

^Different  sample  designs  and  sampling  intensities  were  used  on  the  different  watersheds;  therefore 
the  reliability  of  the  estimates  differ. 


cumstances,  areas  of  uniform  treatment  will  be 
limited  to  several  hundred  acres  in  size,  thus 
making  relationship  B  most  applicable. 

Esthetics. — To  evaluate  esthetic  response  of 
the  Beaver  Creek  treatments  we  have  utilized 
existing  criteria  from  Agriculture  Handbook 
434  (U.S.  Forest  Service  1973)  as  well  as  current 
studies  by  Boster  and  Daniel  (1972). 

As  indicated  in  a  previous  section,  the  "en- 
closed," "detail,"  and  "ephemeral"  segments  of 
the  landscape  are  particularly  susceptible  to 
change  by  forest  treatments.  These  changes 
may  be  in  terms  of  "variety,"  "deviations,"  or 
"contrast,"  and  will  result  in  changed  forms, 
colors,  or  textures  of  the  landscape. 

Consider  first  the  effects  of  the  several  land 
treatments  on  "variety."  From  the  standpoint  of 
landscape  quality,  some  variety  is  desirable,  but 
it  is  less  obvious  just  how  much  variety  is 
enough. 

There  is  a  point  at  which  variety,  increasing 
from  zero,  becomes  visually  pleasant  or  accept- 


able. As  variety  continues  to  increase,  it  ap- 
proaches the  point  where  it  is  no  longer  pleasant 
(it  may  decrease  to  zero  again).  In  the  sketches 
from  the  Agriculture  Handbook  (fig.  15),  A  and 
G  tend  to  be  the  least  interesting.  The  inter- 
mediate stages  might  prove  to  be  the  most  en- 
joyable. Although  not  everyone  would  rate  each 
sketch  the  same,  most  would  pick  the  sketches 
in  the  middle  range  as  most  interesting  and  C  or 
E,  where  the  proportions  are  not  equal,  are 
likely  to  be  preferred.  Although  a  precise  pro- 
portion cannot  be  determined,  we  can  approxi- 
mate the  range  on  the  variety  scale  at  which 
visual  acceptability  is  reached. 

Figure  15  provides  a  basis  for  evaluating 
changes  in  variety  brought  about  by  several  of 
the  Beaver  Creek  treatments.  We  might  infer 
that  neither  the  untreated  nor  the  clearcut  areas 
provide  desirable  variety.  The  strip  shelter- 
wood  cut  on  Watershed  14  would  have  the  high- 
est rank  of  the  treatments.  The  severe  thinning 
and  severe  shelterwood  treatments  on  Water- 
sheds 17  and  16  would  have  intermediate  vari- 
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Table  7- ""Summary  of  effects  of  treatments  on  deer  and  elk  pellet  group 
densities  in  the  pine  type  on  Beaver  Creek^ 


Watershed,  by  treatment, 
and  year 


Pe 11 et  g  roup 
dens  i  ty , 


Difference  between  actual  and  predicted 
pellet  group  densities  by  years 
following  treatment^ 


Eval uation^ 


treatment  completed 

unt  reated 
cond  i  t  i  ons 

1st 

2d 

3d 

kth  5th 

Aver- 
age 

- 

-  -  - 

Pellet  group 

per 

acre 

per  month  -  - 

Watershed  12-- 

Dee  r 

1.2 

1 

.3 

1.6 

0  1.6 

1  . 1 

Meaningful 

100%  clearcut 

1967 

Elk 

1  . 0 

-.6 

- 

.6 

-.2 

(-)  .2 

-.2 

Not  meaningful 

Watershed  9"" 

Dee  r 

.5 

-.1 

.9 

.2 

".  1 

.2 

Not  meaningful 

32%  clearcut  in 

Elk 

1.2 

-1.0 

- 

.6 

-.8 

1.2 

-.3 

Not  meaningful 

uniform  strips 

1968 

Watershed  1 1 -- 

Deer 

1     0  ■ 

-.6 

.  6 

■1 .0 

-.8  -1.1 

-.8 

Not  meaningful 

60%  ut i 1 i zat ion  of 

Elk 

11.8 

-8.6 

-10 

.i*  ■ 

■8.7 

-10.6  -10.8 

-9.8 

Meaningful 

\JG  iciiiiiai  Ljia^d 

1  Q(i7 
1  / 

Watershed  17"" 

Deer 

.8 

.2 

3 

.6 

1.9 

1.9 

Mean  i  ngf ul 

75%  thin 

1969 

Elk 

1.0 

.3 

6 

.3 

3.6 

Meani  ngf ul 

Watershed  l'*-- 

Deer 

.6 

1.5 

1.5 

? 

33%  pine  clearcut  in 

Elk 

.if 

.6 

.6 

? 

irregular  strips, 

thinning  between 

strips;  50%  basal 

area  removed  overall 

1970 

Adapted  from:     Neff,  Don  J.     Effect  of  watershed  treatment  on  deer  and  elk  use.     Job  Progr.  Rep. 

Proj.  W78-16,  Work  Plan  h.  Job  5.  1973- 
^Predicted  production  is  an  estimate  of  production  that  would  be  expected  assuming  the  watershed 

had  not  been  treated. 

^Evaluation  of  the  treatment  effect  based  on  the  judgment  of  the  Arizona  Game  and  Fish  Department 

Biolog  i  St . 
''Trace. 


Figure  1 4.— Relationships  of  deer  use  to  timber  basal  area:  A,  Areas  of  <  40  acres  surrounded  by  forest 
cover  of  3=  1 00  ft.2/acre  of  basal  area;  B,  Areas  of  40  to  500  acres  surrounded  by  forest  cover  of  >  1 00 
ft.2/acre  of  basal  area;  C,  Areas  of  uniform  timber  basal  area. 
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ety  rankings.  If  we  related  the  sketches  to  basal 
area,  it  would  appear  that  the  intermediate 
basal  area,  50  to  80  ft'/acre,  would  have  the 
highest  variety  ranking,  especially  if  the 
residual  stand  is  arranged  in  groups  of  trees. 

The  Handbook  suggests  that  the  intermediate 
sketches  in  fig.  16,  (C  and  D)  provide  the  more 
desirable  degrees  of  variety.  There  are  definite 
analogies  in  the  uniform  stripcut  treatment  on 
Watershed  9  (Sketch  A)  and  the  irregular  shel- 
terwood  stripcut  on  Watershed  14  (Sketch  D). 
The  Watershed  14  treatment  scores  well  in  var- 
iety except  where  the  strips  appear  too  repeti- 
tive. The  uniform  stripcut  on  Watershed  9  is 
much  less  desirable. 

Treatments  that  create  unnatural-appearing 
deviations  from  the  surrounding  landscape  are 
a  distinct  esthetic  detraction.  Watershed  12, 
which  was  clearcut,  is  in  this  category  since  it  is 
essentially  a  single  large  opening  covered  with 
parallel  windrows  of  slash.  An  irregular  pattern 
of  forest  openings  which  borrow  form,  line, 
color  and  texture  from  the  surrounding  land- 
scape would  be  preferable. 

The  ephemeral  quality  of  the  landscape  can 
be  enhanced  by  increasing  the  proportion  of 
Gambel  oaks  and  aspen  in  forest  stands.  Gambel 
oak  has  been  retained  on  Watersheds  14, 17,  and 
16  primarily  for  its  wildlife  values,  with  the 
added  result  that  fall  colors  are  substantially 
enhanced. 

Esthetics  are  being  measured  on  ponderosa 
pine  areas  of  the  Coconino  National  Forest  and 
elsewhere  by  means  of  an  extension  of  the  TSD 
(theory  of  signal  detection)  to  evaluate  viewer 
preference  for  alternative  land  treatments.  The 
TSD  model  appears  to  have  considerable 
promise  for  this  type  of  evaluation  (Boster  and 
Daniel  1972;  Boster  1973;  Daniel  et  al.  1973). 

TSD  methodology  was  designed  principally  to 
distinguish  between  the  two  components  of  per- 
ceptual judgment:  actual  perception  and  the  ob- 
servers' standards.  As  noted  in  Daniel  et  al. 
(1973),  it  was  further  developed  in  response  to  a 
need  to  measure  image  quality  when  complex 
photographic  displays  are  subjected  to  infor- 
mation retrieval  or  enhancement  processes. 
The  Forest  Service  is  attempting  to  use  TSD  to 
answer  two  questions  related  to  forest 
management:  (1)  can  observers  discriminate 
reliably  among  various  vegetative  treatments, 
and  (2)  do  observers  make  differential  esthetic 
responses  to  the  various  treatments?  Ulti- 
mately the  TSD  approach  permits  the  quantifi- 
cation of  the  esthetic  judgments  on  a  cardinal 
(or  interval)  scale  rather  than  on  merely  an 
ordinal  (ranking)  scale. 

During  its  developmental  stage  this  evalua- 
tion technique  was  tested  extensively  with 
selected  experimental  watershed  treatments. 


A 


Figure  15. — Degrees  of  variety  in  a  forest  landscape 
(from  U.S.  Forest  Service  1973). 


After  sufficient  validity  and  reliability  testing, 
a  different  set  of  areas  was  selected  represent- 
ing a  wide  array  of  feasible  management  op- 
tions. Numerous  definable  groups  of  both  forest 
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Figure  16. — Degrees  of  variety  in  the  forest  land- 
scape along  a  road  or  other  cut  strip  (from  U.S. 
Forest  Service  1973). 

users  and  nonusers  are  now  being  questioned 
about  their  preferences  for  these  alternative 
management  practices: 


Management  practice 
Uncut^ 

Irregular  stripcut 
Heavy  thin 

Moderate  thin 

Recent  management 
Intensive  management 


Area  Description 

Woods  Canyon  Decision 
Area 

Watershed  14 
Watershed  17,  30  ftV 
acre  basal  area 
Wild  Bill  Pasture,  60  ft^/ 
acre  basal  area 

'69  NFS  "Horse  Lake" 
sale  area 

72  NFS  "Fred  Hought 
Ridge"  sale  area 
(a  silvicultural 
prescription) 


One  of  the  most  critical  esthetic  considera- 
tions concerns  how  cutting  changes  scenic  qual- 
ity. There  is  general  agreement  that  scenic 
quality  drops  markedly  immediately  following 
harvest,  then  increases  as  vegetation  recovers 
and  logging  scars  heal.  The  main  question  is. 
How  much  recovery  can  we  expect?  Related 
questions  concern  whether  the  pretreatment 
scenic  quality  (A,  fig.  17)  will  be  reached  or 
exceeded  (curves  a'  and  a,  fig.  17)  and  how  much 
time  will  be  required  (x  years  in  fig.  17).  Ans- 


Judged 
scenic 
quality 


Figure  1 7. — Esthetic  recovery  of  an  area  harvested  in 
yearly.  See  text  for  explanation  of  symbols. 

wers  to  these  questions  are  being  sought  by  Bos- 
ter  and  Daniel  for  several  recently  treated 
areas.  The  "recovery"  of  Watersheds  9, 14,  and 
17,  plus  that  of  a  recently  logged  area  that  was 
lightly  cut  (Milk  Ranch  Point),  have  been  fol- 
lowed for  several  years.  These  answers  can  be 
expected  to  go  a  long  way  toward  filling  an  im- 
portant void  in  forest  management,  providing 
scientific  justification  for  various  management 
options  as  they  relate  to  scenic  quality. 

Boster  and  Daniel  have  conducted  an  exten- 
sive analysis  to  identify  those  factors  which 
consistently  increase  or  decrease  viewer  satis- 
faction. This  analysis  has  resulted  in  a  rather 

'The  "uncut"  area  that,  although  lightly  cutover,  could 
lend  itself  to  any  of  the  five  subsequent  management  prac- 
tices. 
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small  list  of  significant  factors.  It  has  also  fos- 
tered development  of  a  management  tool  that 
will  enable  landscape  architects  and  others  to 
predict  scenic  quality  solely  from  viewer  judg- 
ments of  landscape  features.  The  factors  used 
must  be  periodically  updated  because  people's 
tastes  and  preferences  change  over  time. 

Landscape  architects  in  Region  3  are  being 
asked  to  scale  representative  slides  with  those 
factors  found  to  best  predict  judged  scenic  qual- 
ity. Based  on  their  scaling,  simple  regression 
equations  will  be  derived  with  scenic  quality  as 
the  dependent  variable  and  a  small  number  of 
identifiable  landscape  features  as  the  depen- 
dent variables.  These  equations  will  be  usable 
by  the  landscape  architects.  The  next  step  will 
be  to  predict  judgments  of  scenic  quality  from 
preharvest  prescriptions  and  sketches  of  pro- 
posed treatments  (or  from  photos  of  existing 
stands  typical  of  the  proposed  harvest  prescrip- 
tion). 

The  relation  between  most  of  the  factors  that 
affect  a  viewer's  satisfaction  and  the  scenic 
value  estimator  statistic  shows  good  agreement 
with  theory.  For  example,  variety  and  diversity 
(as  represented  by  "tree  density"  and  "average 
tree  diameter")  were  found  to  be  important  fac- 
tors. Also,  scenes  lacking  harmony  with  the 
characteristic  landscape  (those  having  large 
amounts  of  down  wood  or  slash,  for  example) 
have  low  preference  ratings.  Very  open  or  very 
dense  forest  stands  likewise  were  less  pre- 
ferred than  stands  with  intermediate  basal  area 
levels.  Figure  21  (page 30)  illustrates  this  latter 
relationship,  with  esthetic  preference  plotted  as 
a  function  of  timber  basal  area.  These  and  other 
preliminary  relationships  are  being  quantified 
for  use  in  a  model  that  will  predict  esthetic  pref- 
erence on  the  basis  of  a  wide  variety  of  perti- 
nent variables. 

Treatment  Costs 

Since  treatment  costs  are  covered  in  detail  in 
another  paper  (Turner  and  Larson  1974),  only  a 
brief  summary  of  costs  will  be  presented  here. 
Turner  and  Larson  have  developed  regression 
models  which  predict  thinning  and  piling  costs 
as  a  function  of  the  degree  of  timber  basal  area 
removed.  Thinning  costs  are  related  to  basal 
area  removed  noncommercially,  while  piling 
costs  are  related  to  total  basal  area  removals 
including  commercial  logging.  The  costs  are 
predicted  in  terms  of  dollars,  and  of  man  and 
equipment  hour  requirements  per  acre,  so  that 
wages  and  rates  can  be  updated  for  any  situa- 
tion. 

Sensitivity  analyses  indicate  that  these  piling 
costs  are  representative  for  all  but  the  most 
extreme  conditions  of  slope  steepness  likely  to 


be  encountered  in  the  Southwest.  If  the  thinning 
operation  involves  removal  of  trees  larger  than 
8  inches  in  diameter,  cost  variability  may  be 
greater  than  that  accounted  for  by  the  thinning 
model. 

By  basing  estimates  on  thinning  and  piling  op- 
erations, which  account  for  about  90  percent  of 
the  total  costs,  we  can  eliminate  variation  due  to 
miscellaneous  lesser  jobs  such  as  travel  which 
varies  between  watersheds. 

A  general  model  applicable  to  both  experi- 
mental and  forest  operations  was  obtained  by 
relating  costs  to  basal  area  removed.  The  cost  to 
remove  noncommercial  trees  was  $0.79  per  ft^ 
of  basal  area,  while  the  piling  cost  was  $0.28  for 
each  ft-  of  commercial  and  noncommercial 
basal  area  removed. 

These  costs  are  based  on  initial  treatments. 
Estimates  of  total  costs  require  projections  to 
include  subsequent  precommercial  thinnings 
and  harvests  throughout  a  complete  rotation. 
Because  treatments  with  high  initial  costs  may 
have  lower  costs  at  future  cuts,  evaluation  of  a 
treatment  alternative  on  the  basis  of  initial 
costs  could  be  very  misleading  unless  future 
product  flows  and  costs  are  accounted  for  and 
discounted  back  to  present.  This  procedure  can 
be  handled  internally  by  a  multiple  use  alloca- 
tion model  which  considers  costs  and  yields  at 
appropriate  intervals  throughout  the  planning 
horizon  (Turner  1974). 

Preliminary  Models  for  Predicting  Resource 
Responses  in  the  Ponderosa  Pine  Type 

A  long-range  aim  of  the  Forest  Service  re- 
search program  in  Arizona  is  to  provide  an 
array  of  analytical  procedures  (models)  for 
predicting  the  following  responses  to  manage- 
ment activities: 

1.  Terrestrial/biotic  responses  (timber,  range, 
and  wildlife). 

2.  Hydrologic  responses  (water  yield,  quality, 
and  timing,  erosion,  and  sedimentation). 

3.  Stream  and  lake/biotic  responses  (fish  habi- 
tat and  production). 

4.  Natural  beauty  and  recreation  responses 
(esthetic  value  and  recreation  potential). 

5.  Fire  responses  (fire  size  and  rate  of  spread). 

Response  data  from  such  models  will  provide 
input  to  allocation  and  scheduling  models, 
which  can  be  used  to  allocate  uses  and  activities 
to  areas  of  land  in  a  planning  unit  and  to 
schedule  activities  on  all  planning  units  in  a 
forest  or  region. 

Physical  models  are  already  available  for 
predicting  some  of  these  responses  on  pon- 
derosa pine  forest  lands.  A  continuing  effort 
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will  be  made  to  determine  the  suitability  of 
these  models  and  to  utilize  them  to  the  extent 
possible  in  building  an  overall  system  of  mutu- 
ally compatible  models  suitable  for  use  in  plan- 
ning and  managing  forest  and  related  lands  in 
the  Southwest.  This  overall  system  of  models 
wUl  eventually  be  expected  to  produce  all  of  the 
response  information  listed  above. 

Of  the  existing  water  yield  models  available, 
the  Stanford  Watershed  Model  (Crawford  and 
Linsley  1966)  is  not  suited  for  application  on 
small  watersheds.  Because  the  USDA  model 
(Holtan  and  Lopez  1971)  is  not  capable  of  ac- 
counting for  snowmelt,  it  is  not  applicable  to 
areas  such  as  the  ponderosa  pine  type  where  a 
majority  of  runoff  is  derived  from  a  melting 
snowpack.  Water  Yield  I  (U.S.  Forest  Service 
1972)  may  be  applicable  to  the  Beaver  Creek 
area  but  is  not  yet  available  for  testing.  Leaf  and 
Brink's  (1973)  hydrologic  model  will  be  tested 
once  it  has  been  modified  to  function  in  areas 
with  shallow  or  intermittent  snowpacks.  Test- 
ing of  the  BURP  model  (U.S.  Forest  Service 
1968)  is  also  being  planned.  Such  testing  work 
has  been  delayed,  however,  pending  completion 
of  an  overall  problem  analysis  concerning  the 
development  of  models  needed  to  provide 
management  information  for  the  Southwest 
forest  ecosystems. 

Meanwhile,  preliminary  models  can  be  used 
to  provide  some  initial  arrays  of  response  data. 
Two  water  yield  models  have  undergone  pre- 
liminary testing  on  Beaver  Creek:  Baker  and 
Kovner's  regression  model,  and  Rogers'  system 
theoretic  model. 

Baker-Kovner  Streamflow  Regression  Model 

First  attempts  at  modeling  streamflow  on 
Beaver  Creek  have  used  a  regression  approach. 
A  multiple  regression  to  predict  annual 
streamflow  was  developed  from  148  observa- 
tions from  12  watersheds  in  the  pine  type.  Vari- 
ables initially  considered  included  precipita- 
tion, insolation,  soil,  geology,  and  timber  den- 
sity. 

The  objective  was  to  formulate  a  model  that 
was  in  accord  with  physical  concepts  influenc- 
ing the  hydrologic  regime  on  the  watersheds. 
Only  parameters  which  logically  influenced  the 
water  balance  were  considered.  Relationships 
between  parameters  were  also  analyzed  to  de- 
termine if  they  were  independent,  linear,  or 
curvilinear. 

The  initial  model  was  a  linear  multiple  re- 
gression which  contained  six  variables.  This 
equation  was  cumbersome,  and  examination 
showed  that  not  all  variables  were  independent. 
It  also  became  apparent  that  the  relationship 
was  nonlinear.  A  following  modification  of  the 


regression  included  winter  precipitation,  po- 
tential insolation,  and  timber  density.  It  too  was 
nonlinear. 

Since  93  percent  of  the  runoff  on  Beaver 
Creek  occurs  in  the  winter,  winter  precipitation 
was  found  to  correlate  best  with  annual 
streamflow.  The  potential  insolation  variable 
being  used  is  a  percentage  of  the  potential  inso- 
lation falling  on  a  surface  normal  to  incoming 
radiation  (Fons  et  al.  1960).  At  present  this  vari- 
able has  been  indexed  by  using  the  present 
normal  insolation  occurring  at  1200  hours  on 
February  23.  This  particular  date  was  selected 
because  much  of  the  snowmelt  on  Beaver  Creek 
occurs  around  this  time.  The  insolation  variable 
may  be  improved  in  the  future  by  using  average 
annual  or  winter  insolation  values  for  the  vari- 
ous watersheds. 

Although  the  relationships  developed  with  the 
basic  data  from  the  Beaver  Creek  watersheds 
between  streamflow,  winter  precipitation,  and 
insolation  are  linear,  the  expected  relationship 
between  vegetation  and  streamflow  is  S-shaped 
because  of  known  upper  and  lower  limits. 
Timber  density  parameters  considered  were 
number  of  trees,  basal  area,  and  cubic-foot  vol- 
ume per  acre.  Of  these,  basal  area  and  cubic- 
foot  volume  appear  to  be  most  useful. 

There  is  obviously  interaction  between  pre- 
cipitation, insolation,  and  timber  density.  Two 
interaction  terms  used  in  the  regression  are  be- 
tween winter  precipitation-basal  area  and  be- 
tween winter  precipitation-insolation.  The  form 
of  the  regression  which  has  been  developed  is: 

Y  =  -5.72  +  0.83Xi  +  0.42X2  -  0.24X2X1°'' 

-0.007Xi'  [1  —  exp  —  (X3/45)]3 
R'  =  0.69 

Y  =  Annual  streamflow  in  inches 
Xi  =  Winter  precipitation  in  inches 
X2  =  Insolation  as  a  decimal  fraction 
X3  =  Timber  basal  area  in  ftVacre 

The  observational  data  for  such  a  model  is 
limited  because  there  are  only  12  watersheds. 
On  Beaver  Creek,  winter  precipitation  has 
ranged  from  8  to  28  inches,  insolation  from  0.66 
to  0.74,  timber  volume  from  0  to  2,500  ft^/acre, 
and  timber  basal  area  from  0  to  125  ft^/acre.  It  is 
difficult  to  accurately  control  the  estimates  of 
eight  coefficients  in  linear  and  nonlinear  form 
over  the  full  range  of  combinations  with  a  small 
amount  of  data.  Therefore  the  parameters  are 
not  necessarily  functionally  related  to  the 
physical  processes  which  control  the  water 
balance. 
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In  figure  18  the  model  provides  graphical  es- 
timates of  streamflow  for  three  winter  pre- 
cipitation levels  and  a  range  of  basal  area  levels. 
To  illustrate  use  of  the  curves,  we  would  esti- 
mate that  streamflow  will  be  increased  0.6  inch 
if  basal  area  is  reduced  uniformly  from  120  to  60 
ft-/acre  in  a  situation  where  winter  precipita- 
tion is  15  inches.  The  model  predicts  4.3  inches 
of  streamflow  as  an  average  for  watersheds  on 
Beaver  Creek  with  16  inches  of  winter  precipi- 
tation and  basal  area  of  116  ft-/acre.  Actual 
streamflow  from  these  same  watersheds  is  5.3 
inches.  For  estimating  response  to  strip  or  patch 
cutting,  it  is  probably  best  to  consider  the  clear- 
cut  openings  to  be  zero  basal  area,  and  use 
actual  basal  area  for  remaining  timber  stands. 

Since  on  the  average,  only  5  years  of  post- 
treatment  data  are  available  from  Beaver 
Creek,  predicted  responses  should  only  be  used 
to  estimate  average  response  for  a  limited 
period  following  treatment.  This  average  re- 
sponse will  decrease  with  time  after  treatment, 
and  is  estimated  to  be  of  no  practical  impor- 
tance after  20  years. 

Over  a  14-year  period  on  three  watersheds, 
the  model  is  predicting  about  20  percent  below 
derived  streamflow  (table  8).  When  all  observa- 
tions from  the  12  pine  watersheds  are  con- 
sidered 69  percent  of  the  variation  in  streamflow 
is  accounted  for  by  the  regression,  with  the  ma- 
jority being  accounted  for  by  winter  precipita- 
tion. The  regression  model  generally  under- 
estimates yields,  but  is  relatively  consistent 
even  in  an  unusually  high  runoff  year  such  as 
1973. 

The  model  was  also  tested  on  data  from  the 
East  and  West  Fork  ponderosa  pine  watersheds 
on  Castle  Creek  in  east-central  Arizona  (table  9). 
East  Fork,  which  is  used  as  the  control,  has  a 
timber  basal  area  of  122  ft'/acre.  Its  insolation  is 
estimated  to  be  0.61.  West  Fork,  where  timber 
was  commercially  harvested  by  prescription, 
had  pre-  and  posttreatment  timber  basal  areas 
of  135  and  63  ft'/acre,  respectively.  The  insola- 
tion on  West  Fork  is  estimated  to  be  0.72  (Rich 
1972). 

The  16-year  average  actual  streamflow  on 
East  Fork  is  2.8  inches,  compared  with  a  pre- 
dicted average  of  2.7  inches  (table  9). 

The  pretreatment  11-year  average  actual 
streamflow  on  West  Fork  is  2.2  inches  and  the 
predicted  mean  is  2.9  inches;  the  posttreatment 
5-year  average  of  actual  streamflow  is  2.0  in- 
ches and  the  predicted  mean  is  1.9  inches. 

The  regression  model  predicted  average  pre- 
treatment streamflow  for  West  Fork  of  Castle 
Creek  within  32  percent  and  posttreatment 
streamflow  within  6  percent.  The  long-term 
mean  —  including  pre-  and  posttreatment 
periods  —  was  estimated  within  21  percent  of 


the  actual.  The  model  estimated  a  streamflow 
response  of  about  0.3  inch,  whereas  the  "actual" 
response  determined  by  averaging  individual 
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Figure  18. — Predictions  of  annual  streamflow  by  re- 
gression analysis  based  on  winter  precipitation, 
insolation,  and  timber  basal  area. 


Table  8. --Compar i son  of  observed  and  predicted 
streamflow  (inches)  from  three  water- 
sheds on  Beaver  Creek 
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Table  9 . --Compa r i son  of  observed  and  predicted 
streamflow  (inches)  from  East  Fork 
and  West  Fork  watersheds  on  Castle 
Creek  in  eastern  Arizona 
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yearly  responses  was  0.5  inch.  The  regression 
estimate  for  the  combined  pre-  and  posttreat- 
ment  periods  on  the  EastFork  of  Castle  Creek 
was  within  4  percent  of  the  actual. 

This  limited  test  indicates  that  average 
streamflow  can  generally  be  estimated  within  20 
or  30  percent  of  the  actual  in  the  ponderosa  pine 
type  with  volcanic  soils  and  climate  similar  to 
those  on  Beaver  Creek.  Error  in  streamflow  es- 


timates for  individual  years  may  be  larger,  and 
the  model  cannot  be  used  on  the  more  perme- 
able sedimentary  and  alluvial  soils  without 
some  adjustment. 

Rogers'  System  Theoretic  Hydrologic  Model 

Rogers  (1973)  has  developed  a  model  of  the 
hydrologic  behavior  of  a  managed  natural 
ecosystem,  following  concepts  published  by 
Wymore  (1967,  1972).  The  model  assumes  that 
an  ecosystem  may  be  divided  into  relatively 
homogeneous  land  units  and  channel  reaches.  It 
requires  as  inputs  inventory  data  describing 
physiographic  and  vegetative  conditions  on 
each  unit  and  reach,  along  with  information  on 
climate  and  treatments.  It  operates  on  a  time 
interval  varying  from  5  minutes  to  12  hours. 

The  water  and  energy  balance  are  computed 
on  each  unit  of  the  ecosystem.  Water  yields  are 
routed  from  unit  to  unit  and  into  the  channel 
system.  Major  processes  modeled  include  those 
affecting  canopy  interception,  snowpack  water 
and  energy  balance,  mulch  layer  water  balance, 
surface  water  balance  and  overland  flow  rout- 
ing, evapotranspiration,  soil  water  balance,  in- 
terflow routing,  and  channel  routing. 

Rogers  (1973)  includes  a  computer  program 
for  implementing  the  model,  together  with  in- 
structions for  its  use.  He  also  gives  a  procedure 
for  developing  the  required  input  data. 

The  model  was  independently  tested  on 
paired  Watersheds  12  and  13.  Twelve  years  of 
data  were  available.  Watershed  12  had  been 
clearcut  after  the  seventh  year  while 
Watershed  13  was  a  control.  The  predicted 
mean  (5.27  inches)  and  standard  deviation  (4.85) 
of  annual  flows  for  the  clearcut  watershed  was 
approximately  the  same  as  the  observed  (5.42 
and  4.27,  respectively,  table  10).  The  predicted 
mean  (1.96)  and  standard  deviation  (1.91)  for  the 
untreated  watershed  were  significantly  less 
than  the  observed,  however  (3.10  and  3.08,  re- 
spectively). 

Comparisons  of  actual  and  predicted  hydro- 
graphs  (fig.  19)  indicate  that,  in  general,  the  pre- 
dicted daily  flows  follow  the  same  general  pat- 
tern as  the  observed  flows.  During  the  1960, 
1961,  1963,  1964,  1965,  and  1966  calendar  years 
the  model  generally  predicted  the  start  of 
snowmelt  period  and  time  of  peak  flows  to 
within  several  days.  However,  the  volumes  of 
peak  flows  were  generally  outside  the  error 
limits.  In  1959  the  December  melt  period  was 
not  predicted,  and  in  1962  the  early  melt  period 
was  missed  although  later  ones  were  predicted. 
In  1969  the  time  of  major  peaks  was  predicted 
although  the  smaller  ones  were  missed,  and  in 
1971  most  of  the  later  peaks  were  predicted  al- 
though the  first  one  was  missed.  For  the  150-  to 
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Table  1 0. --Compar i son  of  observed  and  predicted 
streamflow  (inches)  from  two  water- 
sheds on  Beaver  Creek  (Rogers'  system 
theoretic  hydrologic  model) 
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200-year  Labor  Day  Storm  of  1970,*  the  model 
overestimated  the  actual  flows. 

It  is  difficult  to  evaluate  the  model's  perform- 
ance with  these  limited  test  results. 

Primary  sources  of  error  appeared  to  be  in  the 
climatic  inputs,  some  of  which  had  to  be  esti- 
mated, and  in  modeling  the  various  effects  of 
treatments  on  soil  and  vegetation  conditions. 
More  comprehensive  testing  is  needed. 

Timber  Yield  Simulation  Model 

The  first  attempt  at  timber  yield  modeling  on 
Beaver  Creek  was  a  regression  model  de- 
veloped from  the  pine  watershed  data.  The  eight 
independent  variables  tested  included  annual 
precipitation,  insolation,  soil  depth,  cubic-foot 
volume,  winter  precipitation,  elevation,  percent 
cinders,  and  trees  per  acre.  An  equation  using 
the  first  four  independent  variables  for  estimat- 
ing cubic-foot  volume  growth  accounted  for  72 
percent  of  the  variation.  This  predictive  ability 
may  be  improved  by  stratifying  the  watersheds 
by  density  and  site  classes  instead  of  using  a 

"Baker,  MalchusB.,  Jr.,  Harry  E.  Brown,  and  Norman  E. 
Champagne,  Jr.  Hydrologic  performance  of  the  Beaver 
Creek  watersheds  during  a  100-year  storm.  Paper  pre- 
sented at  Amer.  Geophys.  Union  meeting,  San  Francisco, 
Calif,  19  p.  Dec.  1971. 
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Figure  19.— Test  of  the  Rogers'  system  theoretic  hydrologic  model  (Watershed  12,  1965). 
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watershed  average.  Regression  models  will 
never  simulate  actual  stand  development,  how- 
ever, because  they  are  unable  to  accurately  dup- 
licate all  of  the  functions  representing  regener- 
ation, growth,  and  mortality. 

Due  to  the  limitations  of  regression  models,  a 
simulation  model  is  now  being  developed  which 
employs  functions  for  seed  dispersal,  germina- 
tion, seedling  survival,  tree  growth,  competi- 
tion, and  mortality This  model  is  primarily  for 
unevenaged  stands,  but  can  be  used  to  estimate 
future  stocking  of  evenaged  stands  as  well.  The 
model  simulates  growth  of  trees  through 
twenty-two  2-inch  diameter  classes  including  a 
zero  diameter  class  for  trees  less  than  1  inch  in 
diameter,  and  an  establishment  class  for  regen- 
eration. 

The  model  is  initialized  for  each  stand  by 
reading  site  index,  stand  density  level  (basal 
area)  of  residual  stand  following  harvests,  time 
(years)  of  harvests,  type  of  harvests,  and  the 
initial  stocking  table. 

Numerous  functions  control  regeneration  in 
the  model.  Cone  production  is  a  function  of 
stand  density  and  number  of  trees  per  acre  in 
cone-producing  size  classes.  Seed  production  is 
a  function  of  the  number  of  cones  and  the 
number  of  seeds  per  cone.  Seedling  germination 
and  survival  depend  on  total  crown  cover,  seed 
viability,  and  expected  losses  due  to  rodents, 
birds,  and  other  factors. 

Growth  rate  is  controlled  by  site  index,  stand 
density,  and  tree  size  class.  A  modified  stand 
table  projection  is  used  to  move  a  proportion  of 
the  trees  into  the  next  higher  diameter  class  or 
classes.  The  current  model  uses  a  10-year  time 
step  for  each  growth  period. 

A  parabolic  mortality  function  is  employed  to 
reflect  the  higher  mortality  rate  in  small  and 
large  diameter  classes  and  the  relatively  low 
rate  in  the  intermediate  classes. 

A  height-diameter  curve  is  generated  for 
each  simulation  run.  This  curve  is  a  function  of 
site  index  and  average  growth  rate.  The  gener- 
ated curve  is  used  in  volume  calculations  which 
reflect  the  differences  of  tree  form  due  to  site 
quality. 

Simulation  runs  have  been  made  to  determine 
projected  timber  yields  for  the  two  dominant 
soil  types  on  Beaver  Creek,  over  a  range  of 
stand  density  levels  (fig.  20).  The  board-foot 
volumes  projected  by  the  model  are  consistent 
with  results  from  other  studies  (Schubert  1974; 
Clary  et  all")  in  that  yield  potential  is  maximum 

^Larson,  Frederic  R.  Simulating  growth  and  manage- 
ment of  all-aged  southwestern  ponderosa  pine  stands. 
(Ph.D.  Diss,  to  be  completed  in  1975.) 
^O'Clary,  Warren  P.,  William  H.  Kruse,  and  Frederic  R. 
Larson.  Cattle  grazing  and  wood  production  under  differ- 
ent ponderosa  pine  stand  densities.  Unpubl.  ms.  Rocky 
Mt.  For.  and  Range  Exp.  Stn.,  Flagstaff,  Ariz. 


Figure  20. — Sawtimber  yield  as  a  function  of  basal 
area,  for  two  soil  types  on  Beaver  Creek.  Larson's 
model  (see  text  footnote  9)  was  used  to  simulate 
total  yield  of  intensively  managed  stands  over  a 
1 40-year  period;  this  amount  was  then  averaged  to 
obtain  the  estimates  of  annual  yield. 

at  or  near  the  60-ft^/acre  density  level.  The  pro- 
jections by  Schubert  (1974)  indicate  higher 
production  is  possible  on  better  sites,  but  the 
general  shape  of  the  curve  is  similar  to  those 
generated  by  this  model. 

The  projections  by  Schubert  (1974)  and  Clary 
et  al.'"  indicate  that  annual  cubic-foot  volume 
growth  increases  with  stand  density  until  the 
60-ft2  level  is  reached,  and  then  levels  out 
through  the  balance  of  the  density  levels  tested. 
Larson's  cubic-foot  volume  projections  differ 
from  these  in  that  production  falls  off  as  basal 
area  increases  over  80  ft^/acre.  The  differences 
in  these  projections  may  be  due  to  initial  struc- 
ture. Schubert's  projections  are  based  on  even- 
aged  stands  carried  through  maturity.  Clary's 
projections  are  based  on  average  annual  growth 
rate  of  pole-sized  stands.  Beaver  Creek  stands, 
on  the  other  hand,  are  typified  by  an  abundance 
of  saplings  and  small  poles,  a  shortage  of  small 
sawtimber,  and  an  abundance  of  overmature 
sawtimber. 

A  timber  yield  model  developed  by  Myers 
(1973)  has  also  been  tested.  This  model,  de- 
veloped for  evenaged  stands  in  the  Black  Hills, 
uses  a  regression  equation  to  estimate  future 
diameters  of  evenaged  stands.  Since  regenera- 
tion is  not  a  problem  in  the  Black  Hills,  the 
model  does  not  contain  functions  for  cone  pro- 
duction, seed  production,  germination,  and  sur- 
vival. Growth  rate  is  simply  a  function  of  site, 
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present  diameter,  and  stand  density.  Not  ac- 
counting for  regeneration  problems,  both 
cubic-foot  and  board-foot  volume  projections 
from  Myers'  model  were  about  65  percent 
greater  than  those  for  the  same  site  and  initial 
values  projected  by  Larson's  model. 

Evaluating  Management  Alternatives 

The  previously  described  models  can  now  be 
used  to  provide  limited  assistance  in  land  use 
planning.  With  their  help  the  consequences  of 
alternative  management  practices  can  be  ar- 
rayed in  a  way  that  will  allow  the  manager  to 
judge  which  practices  best  meet  his  objectives. 
Such  arrays  will  eventually  be  prepared  using 
mathematical  programming  methods,  but 
meanwhile  the  procedure  can  be  illustrated 
graphically  with  production  functions.  These 
functions  show  how  the  output  of  a  particular 
product,  such  as  water,  timber,  or  herbage,  re- 
sponds to  increasing  degrees  of  treatments, 
which  in  turn  are  a  function  of  increasing 
amounts  of  variable  inputs,  such  as  labor, 
equipment,  and  materials.  General  concepts  of 
production  functions  as  they  apply  to  multiple 
use  management  were  discussed  by  Lloyd 
(1969),  and  some  preliminary  examples  were 
presented  by  O'Connell  and  Brown  (1972). 

The  functions  use  different  levels  of  timber 
basal  area  as  the  treatment  input.  Thus  to  esti- 
mate response  to  a  uniform  thinning  treatment, 
average  pre-  and  posttreatment  basal  areas  are 
used  to  characterize  the  treatment.  For  strip  or 
patch  cut  treatments  zero  basal  area  is  used  for 
the  cleared  portion  of  the  management  unit,  and 
actual  basal  area  for  the  remaining  portion.  This 
rather  gross  treatment  characterization  admit- 
tedly imposes  certain  limitations  on  the 
functions  —  for  example  they  are  not  sensitive 
to  size,  shape,  or  orientation  of  forest  openings 
or  to  juxtaposition  of  treated  areas.  Nor  do  they 
accommodate  treatment  features  other  than 
forest  density  reductions,  such  as  planting, 
seedbed  preparation,  fencing,  salting,  and  so 
forth.  Nevertheless,  forest  density  manipula- 
tions are  an  important  part  of  most  management 
prescriptions,  and  even  crude  applications  such 
as  these  will  provide  general  management 
guidelines. 

Figure  21  shows  preliminary  single  produc- 
tion functions  for  sawtimber  yield,  herbage 
production,  esthetic  preference,  deer  use, 
streamflow,  peak  stream  discharge,  and  sedi- 
ment as  functions  of  timber  basal  area.  These 
functions  have  been  developed  for  average 
conditions  on  Beaver  Creek.  They  are  of  vary- 
ing reliability  at  this  point  of  development.  The 
best  functions  are  probably  those  for  stream- 
flow  and  herbage,  followed  by  sediment  and 
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Figure  21. — Example  of  a  multiple  use  evaluation  of 
alternative  management  practices  (basal  area 
levels)  designed  to  emphasize  different  uses. 

peak  discharge,  all  of  which  were  developed  by 
regression  analysis.  The  other  curves  were 
handfitted  to  more  limited  data.  Applicability  of 
the  functions  elsewhere  has  not  been  tested  yet 
except  in  the  case  of  the  streamflow  regression 
(see  previous  section). 

Table  11  shows  the  above  information  as  per- 
centage changes.  From  these  data  one  can  de- 
termine at  what  basal  area  level  a  particular 
resource  or  impact  is  maximized,  or  stated 
another  way,  the  level  at  which  to  operate  to 
achieve  a  particular  management  emphasis.  Of 
course,  in  this  context  some  management 
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Table  11.--  Percent  changes  in  products  and  impacts  resulting  from  basal  area  reduction 


Management  emphasis  corresponding  to 
basal  area  reduction  from  120  ft^/acre  to; 


Effect  on 

0 

(Range  and 
St  reamf 1 ow 

40 

(Deer 
use ) 

60 

(Sawt  imber) 

100 
(Esthet  ics) 

(Flood 
ment 

1  20 

and  sed  i - 
cont  rol ) 

BROLL 1 AR  SOILS: 

-  -  -  -  Percent  -  -  -  - 

Sawtimber  (Mbm/acre/yr) 

-100 

67 

31 

0 

Herbage          ( Lbs /acre /yr) 

627 

1  36 

82 

18 

0 

Scenic  gual ity  index 

-1  ^2 

-168 

-183 

9 

0 

Deer  use        (pe'LZet  gvowps/ clots /mo) 

1 92 

317 

1 7 

0 

Streamflow  (vnohes/yv) 

35 

25 

17 

4 

0 

Flood  peaks  (ft^/s/mi^) 

62 

40 

30 

9 

0 

Sediment  (tons/dovs/yr) 

838 

550 

412 

138 

0 

SIESTA-SPONSELLER  SOILS: 

Sawt  i  mbe  r  (Mbm/aare/yr) 

-100 

23 

34 

18 

0 

Herbage  (Ibs/aore/yr) 

991 

182 

1  18 

27 

0 

Scenic  quality  index 

-152 

-168 

-183 

9 

0 

Deer  use        (pellet  groups /acre /mo) 

192 

317 

150 

17 

0 

Streamflow    (inches /yr) 

35 

25 

17 

4 

0 

Flood  peaks  (ft^/s/mi^) 

62 

40 

30 

9 

0 

Sediment        ( tons /acre /yr) 

838 

550 

412 

138 

0 

^Herbage  relationships  similar  to  those  in  figure  21  are  also  available  for  Brolliar  soils, 


emphases  differ  by  only  20  ft^  of  basal  area, 
such  as  flood  and  sediment  control  (120  ft^  of 
basal  area)  and  esthetics  (100  ft^  of  basal  area). 
Such  differences  cannot  be  assumed  to  be  statis- 
tically significant  because  of  the  lack  of  preci- 
sion. Within  these  limitations,  however,  such 
functions  might  still  be  used  to  conceptually 
suggest  optimum  basal  area  levels  for  particu- 
lar management  emphases. 

The  production  functions  also  provide  a  basis 
for  estimating  differences  in  products  resulting 
from  various  levels  of  basal  area.  For  example, 
if  the  basal  area  is  reduced  from  120  to  60  ft^/ 
acre  on  Brolliar  soils,  streamflow  will  be  in- 
creased 0.6  inch,  or  17  percent.  Sawtimber 
production  resulting  from  this  amount  of  basal 
area  reduction  would  increase  an  estimated 
29  fbm/acre/year,  or  67  percent.  Herbage  pro- 
duction would  be  increased  an  estimated  90 
pounds  per  acre,  or  82  percent. 

Physical  product-product  relationships  pro- 
vide another  tool  for  evaluating  management 
responses.  They  show  graphically  how  two 
products  or  environmental  impacts  respond  to 
land  treatments.  If  price  information  is  avail- 
able it  is  possible  to  derive  economic  product- 
product  functions  that  can  be  used  to  determine 
the  most  efficient  combination  of  products. 

The  physical  product-product  curves  are 
formed  by  connecting  points  representing 
paired  outputs  at  a  series  of  basal  area  levels. 


The  curves  in  figures  22  and  23  were  prepared 
from  paired  outputs  from  figure  21. 


*  3 
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0  .05  10  .15 

TIMBER  YIELD  ( 1000  BD.  FT/ A C RE ) 

Figure  22. — Product-product  function  for  streamflow- 
timber.  Tine  curve  shows  relative  yields  of  both 
products  at  different  basal  area  levels.  As  forest 
density  is  reduced  from  120  to  60  ft^,  both  stream- 
flow  and  timber  yields  increase.  Below  60,  stream- 
flow  continues  to  increase  but  timber  yields  decline. 
One  possible  interpretation  is  that  basal  area 
should  be  reduced  to  60  because  both  streamflow 
and  timber  yields  are  increasing  to  that  point. 
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Figure  23. — Product-product  functions  for  esthetics- 
timber  and  deer  use-tinnber.  Esthetics  and  timber 
production  both  increase  as  basal  area  is  reduced 
from  120  to  60.  The  curve  suggests  a  compatible 
operating  level  between  60  and  100.  Deer  use  and 
timber  both  increase  at  basal  areas  down  to  60;  this 
curve  suggests  a  suitable  operating  level  between 
70  and  40. 

The  physical  product-product  curves  can  be 
interpreted  by  observing  the  change  in  the  two 
products  as  basal  area  is  reduced  from  a  start- 
ing level  of  120  ft-/acre.  If  both  products  are 
increasing  as  basal  area  is  reduced,  other  things 
being  equal,  the  manager  would  be  inclined  to 
continue  lowering  the  basal  area.  This  is  the 
complementary  phase  of  the  curve.  When  the 
point  is  reached  where  only  one  product  is 


increasing  and  the  other  is  constant  or  decreas- 
ing, there  is  less  incentive  to  continue  reducing 
basal  area.  It  is  in  this  range  that  the  manager 
would  be  likely  to  establish  the  treatment  level. 
There  is  little  justification  for  operating  in  the 
third  phase  where  both  products  are  decreasing 
(illustrated  in  fig.  22).  The  captions  in  figures  22 
and  23  state  some  of  these  interpretations,  the 
assumptions  being  that  only  two  products  are 
being  considered  at  a  time  and  that  costs  are  not 
limiting. 

Physical  product-product  curves  for  numer- 
ous other  combinations  of  products  or  impacts 
can  be  prepared  and  interpreted  similarly; 
however,  not  all  products  or  impacts  have  the 
same  influence  on  decisionmaking. 

Turner  (1974)  has  shown  in  more  detail  how 
the  relationship  can  be  used  in  the  development 
of  management  prescriptions,  considering 
costs  and  values  of  all  products  and  impacts 
together. 

If  carried  to  the  extreme  there  could  be  dif- 
ferent production  functions  for  all  important 
combinations  of  soil,  vegetation,  climate,  and 
physiography.  Hopefully,  however,  improved 
models  of  ecosystem  behavior  can  be  developed 
that  will  produce  the  same  information  in  a 
more  efficient  manner. 

Meanwhile,  the  regression  models  tested  so 
far  have  quantified  a  limited  number  of  site 
factors  affecting  resource  production  and  en- 
vironmental impacts.  These  are  summarized  as 
follows: 


Products 
and  impacts 


Timber 

Herbage 

Esthetics 

Deer  use 

Streamflow 

Flood  peaks 

Sediment 


Influencing  factors 


Timber  Soil  Winter 


density 

X 
X 
X 
X 
X 
X 
X 


X 
X 


precipi- 
tation 


X 


30-min. 

precipi-  Inso- 
tation  lation 
intensity 


X 
X 


Needed  Research  to  Provide  Better  Manage- 
ment Tools 

Improved  analytical  procedures  and  data  are 

needed  for  predicting  quantity  and  quality  of 
trade  offs  among  major  products,  services,  and 
environmental  values  of  ponderosa  pine  lands 
in  the  Southwest.  Past  research  has  provided 
considerable  information  on  the  effects  of 
selected  management  practices  under  specific 
conditions  on  yields  of  individual  market  prod- 
ucts such  as  water,  timber,  and  forage.  But 
little  information  is  available  on  the  various 


32 


combinations  of  these  products  that  are  attain- 
able from  specific  management  units.  Less  is 
available  on  trade  offs  between  market  prod- 
ucts and  the  many  nonmarket  products  and 
environmental  values  attainable  from  a  given 
area. 

In  this  paper  some  analytical  procedures  have 
been  presented  along  with  accompanying  ar- 
rays of  multiple  use  response  data.  Before  these 
procedures  and  data  can  be  widely  used,  how- 
ever, the  capabilities  of  the  various  predictive 
models  must  be  substantially  enhanced  and 
their  performance  tested  over  a  wide  range  of 
conditions.  Some  particularly  important  needs 
are: 


•  Improved  performance  of  hydrologic  mod- 
els to  allow  better  estimates  of  usable  water 
yield  from  porous  sedimentary  and  alluvial 
soil  types. 

•  Improved  capabilities  for  estimating  hydro- 
logic  responses  of  various  forest  cutting 
patterns  and  size  class  distributions. 

•  Timber  yield  simulation  models  that  are  sen- 
sitive to  soils,  site  index,  soil  water  and  nu- 
trients, and  solar  energy.  Timber  yield  mod- 
els should  be  capable  of  responding  to  a  vari- 
ety of  harvesting  schemes  as  well  as  to  disas- 
ters such  as  insect,  disease,  and  mistletoe 
outbreaks  and  fire.  The  models  will  need  to 
include  tested  functions  for  regeneration, 
growth,  and  mortality. 

•  Functions  to  describe  the  quantity  of  fuels 
for  fire  hazard  prediction. 

•  Improved  capabilities  for  modeling  the  re- 
sponse of  wildlife  habitat  and  populations  to 
alternative  management  practices. 

•  Improved  capabilities  for  estimating  nutri- 
ents and  other  factors  that  affect  site  pro- 
ductivity. 

•  Improved  capabilities  for  predicting  esthet- 
ically  pleasing  landscapes  that  will  be  com- 
patible with  land  treatment  practices  de- 
signed for  other  objectives. 


A  computerized  data  bank  should  be  de- 
veloped to  accommodate  resource  characteris- 
tics data,  climatic  data,  process  or  validation 
data,  and  predictive  output  data  from  alterna- 
tive modeling  systems.  Such  a  data  bank  should 
include  all  available  data  from  the  several  ex- 
perimental watersheds  in  the  southwestern 
ponderosa  pine  type  (Beaver  Creek,  Castle 
Creek,  Thomas  and  Willow  Creeks,  Black  River 


Barometer,  and  Workman  Creek).  Also  included 
should  be  data  from  the  small  homogeneous 
subwatersheds  currently  being  instrumented 
on  Beaver  Creek  and  elsewhere  for  model  test- 
ing and  validation. 

Operational  versions  of  the  above  models 
should  be  developed  in  conjunction  with  Na- 
tional Forest  Systems,  along  with  an  operational 
data  bank  and  inventory  system. 

The  entire  multiple  use  management  plan- 
ming  system  should  be  pilot  tested,  including 
inventory,  data  bank,  response  models,  and 
economic  and  management  models.  Such  pilot 
tests  will  be  necessary  to  improve  the  various 
models  and  procedures,  and  to  show  users  as 
well  as  the  general  public  how  these  manage- 
ment planning  procedures  actually  work. 

Sensitivity  analyses  should  be  performed  on 
the  response  and  management  models  to  de- 
termine which  variables  are  most  important  in 
the  decisionmaking  process,  and  what  gaps  re- 
quire additional  research.  On  the  basis  of  these 
sensitivity  tests,  research  efforts  will  be  reor- 
ganized to  allow  the  most  efficient  and  expediti- 
ous development  of  the  models  into  a  form 
widely  useful  throughout  the  region. 


Conclusions 

In  summary,  four  principal  resource  re- 
sponses are  known  to  result  from  manipulating 
ponderosa  pine  forests  growing  on  volcanic 
soils: 

•  Water  yield  increases  of  1  to  2  inches  per 
year  over  a  5-year  period  have  been  realized 
as  a  result  of  various  intensities  and  patterns 
of  timber  reduction.  On  the  more  productive 
pine  areas,  an  average  increase  of  0.6  inch 
per  year  may  be  realistic,  in  light  of  present 
multiple  use  consideration  and  environ- 
mental constraints. 

•  Preliminary  results  of  a  timber  simulation 
model  suggest  that  harvestable  timber 
volume  growth  may  be  increased  even  after 
a  considerable  reduction  of  timber  basal 
area. 

•  Timber  overstory  reduction  results  in  an 
increase  in  understory  plant  production. 
Herbage  production  can  increase  500 
pounds  per  acre  per  year  on  low  to  moderate- 
ly productive  soils  as  the  result  of  complete 
overstory  removal. 

•  In  general,  deer  and  elk  habitat  are  im- 
proved by  opening  up  the  ponderosa  pine 
canopy  to  enhance  understory  forage 
supply. 
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PESTICIDE  PRECAUTIONARY  STATEMENT 


This  publication  reports  research  involving  pesticides. 
It  does  not  contain  recommendations  for  their  use,  nor 
does  it  imply  that  the  uses  discussed  here  have  been 
registered.  All  uses  of  pesticides  must  be  registered  by 
appropriate  State  and/or  Federal  agencies  before  they 
can  be  recommended. 

CAUTION:  Pesticides  can  be  injurious  to  humans, 
domestic  animals,  desirable  plants,  and  fish  or  other 
wildlife  —  if  they  are  not  handled  or  applied  properly. 
Use  all  pesticides  selectively  and  carefully.  Follow  rec- 
ommended practices  for  the  disposal  of  surplus  pesticides 
and  pesticide  containers. 


